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ABSTRACT 

We present reliable multiwavelength identifications and high-quality photometric redshifts for the 462 X-ray 
sources in the k, 2 Ms Chandra Deep Field-South survey. Source identifications are carried out using deep 
optical-to-radio multiwavelength catalogs, and are then combined to create lists of primary and secondary 
counterparts for the X-ray sources. We identified reliable counterparts for 442 (95.7%) of the X-ray sources, 
with an expected false-match probability of k, 6.2%; we also selected four additional likely counterparts. The 
majority of the other 16 X-ray sources appear to be off-nuclear sources, sources associated with galaxy groups 
and clusters, high-redshift active galactic nuclei (AGNs), or spurious X-ray sources. A likelihood-ratio method 
is used for source matching, which effectively reduces the false-match probability at faint magnitudes com- 
pared to a simple error-circle matching method. We construct a master photometric catalog for the identified 
X-ray sources including up to 42 bands of UV-to-infrared data, and then calculate their photometric redshifts 
(photo-z's). High accuracy in the derived photo-z's is accomplished owing to (1) the up-to-date photometric 
data covering the full spectral energy distributions (SEDs) of the X-ray sources, (2) more accurate photometric 
data as a result of source deblending for « 10% of the sources in the infrared bands and a few percent in the 
optical and near-infrared bands, (3) a set of 265 galaxy, AGN, and galaxy /AGN hybrid templates carefully con- 
structed to best represent all possible SEDs, (4) the Zurich Extragalactic Bayesian Redshift Analyzer (ZEBRA) 
used to derive the photo-z's, which corrects the SED templates to best represent the SEDs of real sources at 
different redshifts and thus improves the photo-z quality. The reliability of the photo-z's is evaluated using the 
subsample of 220 sources with secure spectroscopic redshifts. We achieve an accuracy of | Az|/(l+z) w 1% and 
an outlier [with | Az|/(l+z) > 0.15] fraction of 1.4% for sources with spectroscopic redshifts. We performed 
blind tests to derive a more realistic estimate of the photo-z quality for sources without spectroscopic redshifts. 
We expect there are w 9% outliers for the relatively brighter sources (R < 26), and the outlier fraction will 
increase to « 15-25% for the fainter sources {R > 26). The typical photo-z accuracy is « 6-7%. The outlier 
fraction and photo-z accuracy do not appear to have a redshift dependence (for z ~ 0^). These photo-z's ap- 
pear to be the best obtained so far for faint X-ray sources, and they have been significantly (> 50%) improved 
compared to previous estimates of the photo-z's for the X-ray sources in the « 2 Ms Chandra Deep Field-North 
and « 1 Ms Chandra Deep Field-South. 

Subject headings: cosmology: observations — galaxies: active — galaxies: distances and redshifts — galaxies: 
photometry — X-rays: galaxies 



1. INTRODUCTION 

The Chandra Deep Field-North (CDF-N) and Chandra 
Deep Field-South (CDF-S) are the two deepest Chandra 
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surveys (see iBrandt & Hasingerl l2005l for a review), each 
covering w 440 arcmin^ areas with enormous multiwave- 
length observational investments. Together they have de- 
tected about 1160 X-ray point sources (e.g. , iGiacconi et alJ 
I2002t [Alexander et ani2003t iLuo et alJl2008l hereafter LOS), 
with a sky density of ~ 10000 deg~^ at the limiting flux of 
« 2 X 10"'^ erg cm"^ s"' in the 0.5-2.0 keV band. Most 
(« 75%) of these X-ray sources are active galactic nuclei 
(AGNs), often obscured, at z ~ 0.1^; at faint fluxes, some 
of th e sources are starb urst and normal galaxies at z « 0.1-2 
(e.g.. lBauer etal.ll2004 . 

The CDF-N and CDF-S provide an ideal opportunity to 
study AGN cosmological evolution, AGN physics, the role 
of AGNs in galaxy evolution, and the X-ray properties of 
starburst and normal galaxies, groups and clusters of galax- 
ies, large-scale structures in the distant universe, and Galactic 
stars. However, one prerequisite of such research is to as- 
sociate the X-ray sources with their correct counterparts at 
optical, near-infrared (NIR), infrared (IR), and radio wave- 
lengths, and then to determine their redshifts either spectro- 
scopically or photometrically. The CDF-S was recently ex- 
panded from « 1 Ms to « 2 Ms exposure, and more than 130 
new X-ray sources are detected (LOS). Meanwhile, additional 
deep multiwavelength surveys of the CDF-S have become 
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available, including the VLA radio (e.g., iKellermann et alJ 
I2OO8; Mill er et al.l2008h. Svitzer IR (e.g.. Damen et al. 2010), 
NIR K-hmd (e.g., iGrazian et al.' '2006; Taylor et al. 2009), 
and GALEXW (e.g. Morrissey et al. 2005) surveys; there are 
also a few recent spectros copic surv eys that prov ide updated 
redshift information (e.g.. lVanzella et a l. 2008; Popesso et alJ 
I2OO9I) . It is thus crucial to provide up-to-date X-ray source 
identifications in the CDF-S along with their redshifts to fa- 
cilitate follow-up scientific studies. 

One major uncertainty about X-ray source identifications 
in deep surveys is the contamination from spurious matches, 
e.g., because of the h igh source density in deep optical cat- 
alogs. For example, iBarger et al.l (l2003h were able to as- 
sign optical counterparts to f» 85% of the X-ray sources 
in the « 2 Ms CDF-N. However, they estimated the false- 
match probability to be « 10% at R < 24, and it will rise 
to « 25% at 7? « 24-26. A simple error-circle method was 
used in their work, in which the closest optical source within 
a given matching radius of an X-ray source was selected as 
the counterpart. Recently, a more sophisticated likelihood- 
ratio matching approach was adopted for matching X-ray 
sources in deep Chandra and XMM-Newton surveys to op- 
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tical/NIR sources (e.g., Brusa et al 
l2008t iLaird et all |2009; Aird et af 
ratio technique takes into account the positional errors of both 
the optical and X-ray sources and also the expected magni- 
tude distribution of the counterparts, and thus it can signifi- 
cantly reduce the false-match probability at faint magnitudes. 
Furthermore, performing the source matching at some wave- 
lengths other than optical, e.g., NIR or radio, can also reduce 
the false-match probability, as the X-ray sources may be better 
associated with NIR or radio sources. Combining the multi- 
wavelength X-ray source identifications will effectively lower 
the overall false-match probability. 

As many X-ray sources in deep surveys are optically faint 
(and are therefore challenging targets to obtain spectroscopic 
redshifts), a significant portion of the sources have their red- 
shifts determined photometrically. Photometric redshifts are 
derived by fitting spectral energy distribution (SED) templates 
to the observed broadband photometric data, and their ac- 
curacy is largely controlled by the quality of the data avail- 
able. The previous photometri c redshifts for the X-ray sources 
in the w 2 Ms CDF-N (e.g., Barger et al .ll2003l; 'Cap ak et alJ 
20041) and « 1 Ms CDF-S (e.g.. lWolf et a l. 2004; ZhenJetaD 
2004!) are generally accurate to « 10% with w 15-25% catas- 
trophic redshift failures. The superb and improved multi- 
wavelength coverage in the CDF-S has produced full SED 
sampling for the X-ray sources, from UV to IR, allowing 
determination of photometric redshifts to significantly higher 
accuracy. 

In this paper, we present multiwavelength identifications of 
the CDF-S X-ray sources in the main source catalog of L08. 
Given the counterpart information, we compose a photomet- 
ric catalog for these X-ray sources including up to 42 bands 
of UV-to-IR data. Photometric redshifts are then calculated 
and compared to the latest secure spectroscopic redshifts to 
evaluate their accuracy. The identifications and spectroscopic 
redshifts presented here supersede the counterpart and red- 
shift information provided in L08. In §2 we describe in de- 
tail the likelihood-ratio matching method and our multiwave- 
length identification results. In §3 we present the photometric 
catalog and the derivation of photometric redshifts. The accu- 
racy of the photometric redshifts is appropriately estimated. 
In §4 we discuss the nature of the unidentified X-ray sources. 



and present some future prospects toward source identifica- 
tions. We summarize in §5. 

Throughout this paper, all magnitudes are based upon the 
AB magnitude system (e.g., Oke & Gunn 1983), and we adopt 
the latest cosmology with i/o = 70.5 km s"' Mpc"', JIm = 
0.274, and Q.^ = 0.726 derived f rom the five-year WMAP ob- 
servations dKomatsu et alj|2009l) . 

2. MULTIWAVELENGTH IDENTIFICATIONS OF THE 2 MS CDF-S 
X-RAY SOURCES 

2.1. X-ray and Optical-to-Radio Data 

The X-ray sample we chose is the main Chandra source 
catalog for the « 2 Ms CDF-S, which contains 462 X-ray 
sources (L08). These sources were detected in three X-ray 
bands, 0.5-8.0 keV (full band), 0.5-2.0 keV (soft band), and 
2-8 ke V (hai-d band), with a WAVDETECT false-positive prob- 
ability threshold of 1 x lO"*". The resulting source lists were 
then merged to create the main catalog (see §3.3 of L08 for de- 
tails). The survey covers an area of « 436 arcmin^ and reaches 
on-axis sensitivity limits of « 7.1 x 10"'^, « 1.9 x 10"'^ and 
« 1.3 X 10"'^ erg cm"^ s"' for the full, soft, and hard bands, 
respectively. 

We have visually checked the positions of all the X-ray 
sources in L08 during an X-ray spectral extraction process 
(Bauer et al. 2010, in preparation), using the w 2 Ms CDF- 
S Chandra events file and sometimes also the » 250 ks Ex- 
tende d Chandra Deep Field-South (E-CDF-S; iLehmer et alj 
l2005h events file for sources with large off-axis angles. Al- 
though most X-ray sources in L08 have good positions (no 
positional offsets or offsets < 0.24"), there are five sources 
having positions off by > 1". These are sources 62, 90, 125, 
338, and 456, and we shifted their positions by 1.2", 3.2", 
1.3", 2.8", and 1.5", respectively. Sources 90 and 456 are lo- 
cated at large off-axis angles (11.8' and 9.8'); the other three 
sources are faint (« 21-26 detected counts in the soft band 
with errors of about 8-10 counts), and might be composite 
sources (blended with a weak nearby object). The complex 
point spread functions or morphologies of these five sources 
make it difficult to determine the centroid positions, and their 
positions reported in L08 are poorly defined. '"^ 

We searched for counterparts of the « 2 Ms CDF-S X-ray 
sources in several deep optic al-to-radio catalogs, most of 
which cover the entire CDF-S region. Due to dust absorption 
and redshift effects. X-ray sources in deep surveys usually ap- 
pear brighter at redder wavelengths.'^ Therefore we used the 
following catalogs for identification purposes. 



1. The WFI R-hsnA catalog dGiacconi et all 120021; 
iGiavalisco et al.ll2004h . with a 5 ct limiting AB mag- 
nitude of 27.3. 

2. The GOOBS-S HST version r2.0z z-band catalog 
(IGiavalisco et al.ll2004 "'. with a 5 cr limiting AB mag- 
nitude of 28.2. The GOODS-S survey covers a solid 
angle of « 160 arcmin^ in the center of the CDF-S. 



3. The GEMS HST z-h&M catalog dCaldwell et al.ll2008h . 
with a 5 cr Umiting AB magnitude of 27.3. The GEMS 

Four of these sources (all but source 456) do not have any apparent 
counterparts at their original positions. 

" Counterpart searching has also been performed at several shorter wave- 
lengths (e.g., V and U bands), which yield fewer identifications. 



See http : / / archive ■ s tsci . edu/pub/hlsp/ goods/ catalog_r2/| . 
The r2.0z catalog is based on the v2.0 GOODS-S images, which have signif- 
icantly longer total exposure times in the z band than the v 1 .0 images. 



survey was designed to exclude the GOODS-S region, 
and thus does not c over the entire CDF-S. However, 
ICaldwell et all (1200 8) reduced the GOODS-S vl.O data 
identically to the reduction of GEMS, and the published 
z-band catalog covers the entire CDF-S. 

4. The G OODS-S MUSIC K-hand catalog (iGrazian et alJ 
120061) . with a 90% (w 1.8 a) hmiting AB magnitude of 
23.8. The GOODS-S MUSIC survey covers a solid an- 
gle of w 140 arcmin^ in the center of the CDF-S. The 
V2.0 catalog was used (A. Grazian 2009, private com- 
munication). 

5. The MUSYC /T-band catalog dTavlor et al.ll2009h . with 
a 5 cr limiting AB magnitude of 22.4. This catalog is 
not as deep as the GOODS-S MUSIC K-hmd catalog, 
but it covers the entire CDF-S. 

6. The SIMPLE 3.6 fim catalog dDamen et al.ll2010b . with 
a 5 CT limiting AB magnitude of 23.8.'^ 

7. The y UA 1.4 GHz radio catalog from 'Miller et al] 
(l2008h . We used the 5 a catalog provided by N.A. 
Miller (2009, private communication), with a 5 cr lim- 
iting flux density of w 40 /iJy. Compared to the optical 
positions (WFI R band or GOODS-S z band), the ra- 
dio positions have small systematic offsets in right as- 
cension and declination; since the X-ray positions were 
registered to the optical positions (L08), we manually 
shifted all the radio positions by -0.18" in right ascen- 
sion and 0.28" in declination to remove the systematic 
offsets. 

The basic optical/NIR/IR/radio (hereafter ONIR) catalog 
information is listed in Table [T] Note that the GOODS-S and 
MUSIC catalogs do not cover the entire CDF-S. We also in- 
clude in Table [U the number of sources in each of these cat- 
alogs that are within the w 2 Ms CDF-S region (A^o)- The 
locations of the optical, NIR, and IR bands in a sample SED 
are shown in Figure [T] 

2.2. Matching Method 

Source matching to the X-ray catalog was performed for 
each of the optic al-to-radio catalogs described above. We 
used the likelihood-ratio technique to identify the ONI R coun- 
terparts of the X-ray sources (e.g., Sutherland & S aundersI 
ll992HCiliegi"et a l. 2003; Brusa et al. 2005, 2007). Due to the 
high source densities of the ONIR catalogs and centroid errors 
in the X-ray positions, the simple method of searching for the 
nearest counterpart within a small radius could yield a non- 
negligible number of false matches, usually to faint ONIR 
sources where the background source density is high. The 
likelihood-ratio technique (described in detail below) takes 
into account the positional accuracy of both the ONIR and 
Chandra catalogs and also the expected magnitude distribu- 
tion of the counterparts. It assigns a reliability parameter to 
all the possible counterparts, and it mitigates the effect of false 
matches to faint ONIR sources. 

We used the entire GEMS catalog in tlie source-identification process 
for consistency. However, GEMS counterparts liave a lower priority than 
GOODS-S counterparts when we select primary counterparts for the X-ray 
sources (see j|2.3t 

We also searched for SIMPLE 4.5 fim counterparts, which yielded fewer 
identifications than the 3.6 fim catalog. 



3 

4 2 1 0.4 0.2 0.1 




14.0 14.5 15.0 15.5 



Fig. 1. — The SED of source 312 in LOS along with the selected best-fit 
SED template (a composite starburst/spiral galaxy template), as an example 
to illustrate the photometric bands used for source identification, the coverage 
of the broadband photometric data, and the photometric redshift fitting result. 
EiTor bars for the data points are smaller than or comparable to the symbol 
sizes and thus are not shown. Although this source is not detected in the ra- 
dio band, the VLA 1 .4 GHz radio catalog was also used in the identification 
process (not shown in this plot). The source is not detected in the GALEX 
far-UV band, and the 1 a upper limit is shown. There are 42 filters of photo- 
metric data displayed in the SED. Note that not all of the data points are used 
in the photometric redshift calculation (see j|3.1l for details). The discrepancy 
between some photometric data from different catalogs is probably caused by 
AGN variability. 

For an ONIR candidate with a magnitude m (at a given 
band) at an angular separation r from a given X-ray source,'^ 
the likelihood ratio is defined as the ratio of the probability 
of the ONIR object being the true counterpart of the X-ray 
source and the corresponding probability of th e ONIR ob- 
ject b eing a background, unrelated object (e.g., ICiliegi et al.l 

EqqI, 

LR=^-^. (1) 
n(m) 

where q{m) is the expected magnitude distribution of counter- 
parts, f{r) is the probability distribution function of the angu- 
lar separation, and n{m) is the surface density of background 
objects with magnitude m. We assume that the probability 
distribution of angular separati ons follows a Gaussian distri- 
bution (e.g.. lCiliegi et al.ll2003h : 

^^^>=2i^^p(2S)' 

with the standard deviation a = \J (t1 + <j\, where (Tq and ax 
are the 1 a positional errors of the ONIR and X-ray sources, 
respectively. The 1 a positional errors of the X-ray sources 
were estimated by correlating X-ray sources to bright WFI 
7? -band sources and then examining the positional offsets (see 
§3.3.1 of L08 for details). The 1 a positional errors of the 
ONIR sources are either provided by the catalog or estimated 
from the distributions of positional offsets by matching the 
sources to the GOODS-S z-band catalog. The ONIR posi- 
tional errors and X-ray positional errors are listed in Tables [T] 
andlH respectively. 

The magnitude-dependent surface density of the back- 
ground sources, «(m), was estimated using a sample of ONIR 

" The AB magnitudes for radio sources were converted from the radio flux 
densities (e.g.. lOke &"Guni][T983h , m(AB) = -2.5 log(/^) -48.60, where 
is the flux density in units of erg cm"'^ s"' Hz"' . 
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sources that are at a angular distance between 5" and 30" from 
any of the X-ray sources. We chose sources inside such an- 
nular regions to exclude possible counterparts from the back- 
ground calculation. The expected magnitude distribution of 
counterparts q(m) is not directly observable. To estimate q(m), 
we ran the matching procedure in an iterative way. The match- 
ing results from a previous run were used to calculate q(m) for 
the next run. As a first guess, we selected all ONIR sources 
within 1" of any X-ray source, and we assume these ONIR 
objects are representative of the counterparts for the X-ray 
sources with a small contribution from coincident background 
sources. The number of sources in this sample is A^i, and the 
magnitude distribution of these sources is denoted as total(m), 
which was then background subtracted to derive 



real(w) = total(m) - Trr^Nxnim) , 



(3) 



where Nx is the total number of X-ray sources in the rele- 
vant sample (see Col. 9 of Table [U, and ro = 1". The back- 
ground subtraction generally has a small effect (a few per- 
cent) on the resulting magnitude distribution. See Figure |2] 
for an example of real(m) and n(m) in the GEMS z-band cat- 
alog. Due to the magnitude limits of the ONIR catalogs, we 
are only able to detect a fraction Q of all the true counterparts 
{Q = J"'^ q{m) dm). Thus the expected magnitude distribu- 
tion of counterparts q{m) is derived by normalizing real(m) 
and then multiplying by Q: 



q(m) : 



real(OT) 
X;„,i-eal(m) 



Q 



(4) 



The factor Q was first set to Q = N\ /Nx- The magnitude dis- 
tributions, q{m) and «(m), are smoothed using a boxcar aver- 
age (with width of one magnitude) and then used as input in 
the likelihood-ratio calculation. We note that source blend- 
ing, which is most severe in the 3.6 fim catalog, could affect 
the values of q{in) and n{m) by a tiny fraction and thus in- 
troduce small uncertainties in the matching results. However, 
as we choose counterparts carefully combining all seven cat- 
alogs (see §2.3), this will not affect our final identifications. 

Having computed the values of q(m), fir), and n{m), we 
calculated LR for all the ONIR sources within 5" of any of the 
X-ray sources. We then chose a threshold value for LR (Lth) 
to discriminate between spurious and real counterparts. An 
ONIR source is considered to be a counterpart for an X-ray 
source if its LR value exceeds Lth- In some cases, there will 
be more than one counterpart for an X-ray source within our 
5" limiting radius. If the ratio between the highest and the 
second highest LR values is greater than 5, we keep only the 
counterpart with the highest LR value-"; otherwise, the X-ray 
source is considered to have multiple counterparts in this cat- 
alog. For each counterpart for an X-ray source identified this 
way, we calculate a reliability parameter, Rc, which represents 
the probability of it being the correct identification; 



R.= 



LR 



J:LR + {1-Q) 



(5) 



where the sum is over all the counterparts for this X-ray 
source. 

We believe those counterparts with much smaller LR values are coin- 
cident matches; about 16 such unlikely counterparts were identified and re- 
moved in total, considering all the optical-to-IR catalogs. Because the final 
counterpart for an X-ray source is selected after combining the matching re- 
sults for all the catalogs, the chance of removing a possible true counterpart 
here is negligible. 




Fig. 2. — The expected magnitude distiibutions of the background sources 
(thin dashed histogram; normalized by a factor of Trr^Nx) and counterparts 
(thin solid histogram) for the GEMS z-band catalog in the first trial run. The 
counterparts in the first trail run are GEMS sources within 1" of any X-ray 
source, and they were only used to calculate q(m) and thus LR (see j|2.2t . 
The thick solid and dashed curves are the corresponding smoothed (with a 
boxcar average) distributions that represent real(m) and nr^Nxn{m) (see j|2.2l 
for details). The thick dotted curve shows the final smoothed magnitude dis- 
tribution of the counterparts after three iterations. The magnitude distribution 
of the background sources peaks at fainter magnitude than the distribution of 
counterparts. 

The choice of Lth is not unique, and it depends on two fac- 
tors: first, it should be small enough to avoid missing many 
real identifications; second, it should be large enough to pre- 
vent having many spurious identifications. In t h is paper, we 
followed the discussions in ICiliegi et al.l (l2003h . iBrusa et aP 
(2005, 2007), and Civano et al. (2010, in preparation), and 
we chose a Lth that maximizes the sum of sample complete- 
ness and reliability. The sample completeness parameter (C) 
is defined as the ratio between the sum of the reliability of 
all the sources identified as counterparts and the total number 
of X-ray sources, C = X]^c/A^x- The sample reliability pa- 
rameter (R) is defined as the average of the reliabilities of all 
the sources identified as counterparts. For each run of the 
likelihood-ratio method, we tried several values of Lth and 
chose the one which maximizes the sum of C and R. 

After the first trial run of the likelihood-ratio procedure, 
we calculated real(m) as the magnitude distribution of all the 
sources identified as counterparts in the previous run, and 
we approximated the normalization factor Q to be the com- 
pleteness C of the previous run. For each ONIR catalog, the 
likelihood-ratio procedure was run iteratively until the result- 
ing values of C and R converged. This is a stable process with 
mild perturbations, and the results converge quickly, usually 
in w 2-5 iterations. 

2.3. Matching Results 

A summary of the likelihood-ratio matching parameters and 
results is shown in Table[T] We list the total number of X-ray 
sources that are within the coverage of the ONIR catalog (Nx), 
the number of X-ray sources identified with at least one ONIR 
counterpart (Nd), the number of X-ray sources without iden- 
tifications (A^noId), and the number of X-ray sources identified 
with at least two ONIR counterparts (A^Muiti)- We estimated 
the false-match probability of these identifications using a 
bootstrapping approach. We shifted the X-ray source coordi- 
nates by 5"-53" (with a 2" increment between shifts) in right 
ascension and declination respectively, and recorrelated the 
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shifted sources with the ONIR sources using the hkehhood- 
ratio method. This was done for 100 trials, and the results 
were averaged to estimate the number of false matches. The 
expected number of false matches (A^Paise) and the false-match 
probability for each catalog are shown in Table [T] The false- 
match probabilities range from ss 1-9%. The matching to the 
radio sources yields very secure counterparts, with less than 
one expected false match, largely due to the small background 
radio source density. 

For each catalog, we calculated its counterpart recovery 
rate, defined as the expected number of true counterparts 
(A'lD -A^Faise) divided by the number of X-ray sources (A^x), 
shown in Table [T| Compared to the completeness parame- 
ter C, the recovery rate parameter utilizes the bootstrapping 
results, and is a more appropriate estimate of the catalog 
completeness. The MUSIC NIR and SIMPLE IR catalogs 
have the highest recovery rates, while the VLA radio cata- 
log recovers the fewest counterparts. We also calculated an 
X-ray-to-ONlR fraction (X-O fraction) parameter, defined 
as the expected number of true counterparts (Mo-A^Faise) di- 
vided by the number of ONIR sources (A^o) in the catalog. 
The X-O fraction indicates the degree of association between 
the X-ray and ONIR sources. Only w 1-4% of the opti- 
cal/NlR/lR sources are X-ray sources, whereas w 28% of the 
radio sources are detected by Chandra. A larger X-O fraction 
will lead to more reliable source matching results because of 
the lower background source density. Both the recovery rate 
and X-O fraction parameters are dependent on the catalog 
depth. 

Given the excellent radio-source positional accuracy and 
the secure matching results, we can investigate the accuracy 
of the X-ray source positions using the information for the 
radio counterparts.^' In Figure[3] we show the positional off- 
set between the 94 X-ray sources and their radio counterparts 
versus off-axis angle. The median offset is 0.4", and there 
are clear off-axis angle and source-count dependencies. The 
off-axis angle dependence is mainly due to the degradation of 
the Chandra PSF at large off-axis angles, while the count de- 
pendence is due to the difficulty of finding the centroid of a 
faint X-ray source. Although the results are comparable to the 
LOS study of the positional accuracy based on the matching to 
the WFl /?-band sources (see §3.3.1 of LOS), Figure |3] should 
be considered more representative of the typical X-ray posi- 
tional accuracy, as the radio catalog provides the most-reliable 
X-ray source identification. 

Twenty of the 462 X-ray sources do not have a counterpart 
in any of the ONIR databases, based on the likelihood-ratio 
matching to all seven catalogs. Four (XlDs 72, 165, 234, and 
396) of these 20 sources have a faint ONIR source (a z-band 
source in each case, with z = 25.3-26.7) within w 0.7-1.2" of 
the X-ray position; these ONIR sources were not selected by 
the likelihood-ratio matching due to their small likelihood ra- 
tios {LR < Lth). However, given the relatively small positional 
offsets of these four sources, the lack of any other apparent 
counterparts, and the fact that the likelihood-ratio matching is 
a statistical method and cannot possibly identify all the coun- 
terparts, we manually assign them as the counterparts of the 
relevant X-ray sources (see Col. 13 of Table |2]l. Thus 446 

^' Three of the radio counterparts (those of sources 200, 364, and 368) 
are identified as t he lobe components of multi-component radio sources by 
IKellerman n et al.l ^2008). However, the radio counterparts of sources 364 
and 368 are probably two interacting galaxies based on optical imaging and 
spectrosco pic results (e.g., IVlainier i et al. 2008): these t wo radio sources ar e 
blended in lKellermann et alj EOOSh . but well resolved in lMiller et ai] 120081) . 
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Fig. 3. — Positional offset vs. off-axis angle for sources in the main 
X-ray source catalog of LOS that were matched to the VLA radio sources 
IMiller et al.ll20 08). The open circles represent Chandra sources with < 200 
full-band counts, and the solid dots are Chandra sources with > 200 full-band 
counts. The dotted curve shows the running median of all sources in bins of 
2'. The three outliers on the top right comer of the plot are faint X-ray sources 
with large positional errors (crx ~ 0.5"-0.7"). 

X-ray sources have identifications, including these four addi- 
tional likely cases. Note that the false-match probabilities in 
Table 1 do not apply to these four sources. For the remaining 
16 sources without identifications, there is no ONIR source 
within 1 .3". As the ONIR catalogs we used were all generated 
using red bands, we also visually inspected images in bluer 
bands (e.g., U, B, and V) at the positions of these 16 sources, 
and verified that we did not miss any apparent blue counter- 
parts for them. Eight of these unidentified sources are prob- 
ably related to some other objects (e.g., off-nuclear sources; 
see Col. 1 3 of Table|2]l, and the majority of the other eight are 
likely to be spurious X-ray sources (see discussion in 34.1b . 

Twenty of the sources only have a SIMPLE counterpart; 
however, eight of these SIMPLE counterparts have a corre- 
sponding optical/NlR source nearby (within O.S"), and they 
are likely to be the same source. Therefore we used the more 
accurate optical/NIR positions for these eight counterparts 
(see Col. 13 of Tables- 
Some X-ray sources could have multiple counterparts, ei- 
ther from one ONIR catalog as a result of the likelihood- 
ratio matching (e.g., see the A^Muiti parameter in Table [Til, or 
from different optic al-to-radio catalogs. There are 12 X-ray 
sources with multiple counterparts in at least one of the ONIR 
catalogs; none of them has more than two counterparts in one 
band. For the remaining 434 sources, we consider an X-ray 
source to have multiple counterparts if the angular distance 
between two of the counterparts in different bands exceeds 



3 y (jj + (Tj (so that they appear to be different sources), where 

(Ti and (72 are the positional errors in the relevant ONIR cata- 
logs. There are 60 such sources; five of them have three coun- 
terparts. In all, 72 X-ray sources have multiple counterparts 
(with the other 374 X-ray sources having a unique counter- 
part), and they are flagged in Column 12 of Table |2] 

For each of the 446 X-ray sources with an identification, 
we visually inspected all its counterparts in the appropri- 
ate bands, and chose a primary one. The position of this 
source is the most probable position of the counterpart. We 
chose the primary counterpart from, in order of priority, the 
VLA, GOODS-S, GEMS, MUSIC, WFI, MUSYC, or SIM- 
PLE catalog. This order is chosen based on several related 
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factors: the positional accuracy, angular resolution (to mini- 
mize any blending effects), false-match probability, and cat- 
alog depth. There are a few exceptions: (1) if the optical 
counterpart is clearly extended, i.e., a nearby galaxy, we select 
the optical position instead of the VLA radio position because 
the radio emission could arise from an off-nuclear structure (7 
cases); (2) if the optical counterpart is bright and significant, 
while the radio counterpart is detected at a smaller than 7 a 
significance, we prefer the optical position (4 cases); (3) when 
the counterpart from a lower-priority catalog is much brighter 
and better resolved, or it has a significantly larger reliability 
parameter, we prefer this source as the primary counterpart (5 
cases). The final number of primary counterparts selected in 
each catalog (A^Pii) is listed in Table [T] The mean false-match 
probability for the 442 primary counterparts (excluding the 
the four additional likely counterparts) is then estimated by 
taking the average of the false-match probabilities of individ- 
ual catalogs weighted by the number of primary counterparts 
in each catalog. The mean false-match probability is « 6.2%, 
and it represents the identification reliability for the entire X- 
ray sample (excluding the the four additional likely counter- 
parts). The mean false-match probability also has slight off- 
axis angle and source-count dependencies (due to its relation 
to the X-ray positional accuracy); e.g., it is « 6.0% for sources 
within 4', « 7.0% at > 8', « 9.8% for sources with < 30 full- 
band counts, and « 5.8% for > 30 full-band counts. 

A list of the primary counterparts for the 446 X-ray sources 
(out of 462 total sources) is presented in Table |2] with the 
details of the columns given below; we also list all the sec- 
ondary counterparts (for sources with multiple counterparts) 
in Table [3] with a similar format. 

• Column (1): the X-ray source ID number, which is the 
same as the source number in the main catalog of LOS. 

• Column (2): The 1 a X-ray source positional error ctx- 

• Column (3): the ONIR ID number of the primary coun- 
terpart, taken from the corresponding ONIR catalog 
(see Col. 10). For the 16 X-ray sources that lack a coun- 
terpart, the value is set to "0". 

• Columns (4) and (5): the ONIR J2000 right ascension 
and declination of the counterpart, respectively. This 
position is taken from the corresponding ONIR catalog 
(see Col. 10). X-ray sources that do not have a counter- 
part have this column set to "0". 

• Columns (6): the approximate 1 cr positional error of 
the counterpart, in units of arcseconds, taken from Col- 
umn 6 of Table [U 

• Column (7): the angular distance between the X-ray 
source and its ONIR counterpart, in units of arcseconds, 
set to "0" if there is no counterpart. 

• Column (8): the angular distance divided by the 
quadratic sum of the positional errors (cr = \J + a\), 
i.e., the angular distance in units of a. 

• Column (9): the reliability parameter of this counter- 
part (see Eq. 5). X-ray sources that do not have a coun- 
terpart have a value of "0". 

Although we adopted the same positional error for some catalogs (e.g., 
GOODS-S and GEMS), the deeper catalogs (e.g., GOODS-S) can in principle 
provide slightly better positions for the counteipails. 



• Column (10): the ONIR catalog from which the pri- 
mary counterpart was selected. X-ray sources that do 
not have a counterpart have this column set to "• • •". 

• Column (11): the ONIR magnitude of the primary 
counterpart (see Col. 10 for the band in which the mag- 
nitude was measured). X-ray sources that do not have a 
counterpart have this column set to "0". 

• Column (12): a flag on whether the X-ray source has 
multiple ONIR counterparts, set to "1" for true and "0" 
for false. 

• Column (13): notes on some special sources, including 
the 16 X-ray sources that do not have ONIR identifi- 
cations (cases a-e) and a few other sources (cases f-h). 
(a) NoID: there is no ONIR source within w 1 .3-3.6" of 
the X-ray source position, and it is likely that this X-ray 
source does not have an ONIR counterpart. There are 
eight such sources, (b) Edge: this source lies at the 
edge of the « 2 Ms CDF-S field, with a large PSF size. 
The X-ray photons could come from any of a number 
of galaxies in the vicinity. There is one such source, (c) 
Cluster: the X-ray source may be associated with the 
diffuse emission from a galaxy group or cluster. Galaxy 
groups and clusters in the « 2 Ms CDF-S are identified 
through a wavelet detection of extended X-ray emission 
(Finoguenov et al. 2010, in preparation). There are four 
such sources, which are generally faint and weakly de- 
tected in the X-ray. (d) Star: the source is only detected 
in the soft X-ray band, and is 2.7" away from a bright 
(R= 16) star The X-ray photons are probably largely 
from the star, and the X-ray position was not correctly 
determined by WAVDETECT due to the contamination 
of another possible faint X-ray source nearby. There 
is one such source, (e) ULX: this is an off-nuclear X- 
ray source associated with a nearby galaxy. There are 
t wo such sources. One of them (source 301) is reported 
in lLehmer et al.l (l2006 f). (f) Blend: the primary counter- 
part is blended with an adjacent source, and the position 
of the counterpart is poorly determined (positional un- 
certainty w 1-2"). There is one such source, and the 
primary counterpart is from the SIMPLE catalog (and 
there are no detections of this source at shorter wave- 
lengths), (g) Faint: there is at least one faint ONIR 
source within « 0.7-1.2" of the X-ray source position. 
The likelihood ratio parameter of this ONIR source is 
below the threshold, but we manually assigned it to be 
the counterpart. There are four such sources, (h) Posi- 
tion: the source only has a counterpart from the SIM- 
PLE catalog. However, the SIMPLE counterpart likely 
coincides with a nearby optical/NIR source, and we 
choose this optical/NIR source as the primary counter- 
part because of its more accurate position (see above). 
There are eight such sources. 

The locations of the 446 identified and 16 unidentified 
X-ray sources in the CDF-S region are shown in Figure 21 
along with the coverages of the GOODS-S and MUSIC cata- 
logs, which are the only two catalogs that do not cover the 
entire CDF-S. For the 446 sources with an ONIR identifi- 
cation, the mean and median positional offsets between the 
X-ray sources and their primary counterparts are 0.41" and 
0.32", respectively. The maximum offset is « 1.9". Several 
X-ray sources have large positional offsets (> 1.5") to their 
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Fig. 4. — Positions of the 446 identified (open symbols) and 16 unidenti- 
fied (filled symbols) X-ray sources. The open triangles and circles represent 
sources with spectroscopic and photometric redshifts, respectivel y, an d small 
(large) symbols indicate redshifts smaller (larger) than 1.3 (see j|3.3l for de- 
tails about redshifts). The eight filled circles represent sources that are proba- 
bly related to a star, galaxy groups and clusters, off-nuclear sources, or any of 
a number of galaxies at the edge of the CDF-S field (see Col. 13 of Table[2). 
The six filled squares and two filled stars represent the "NoID" (Col. 13 of 
Table[2) sources that were detected at a WAVDETECT false-positive threshold 
of 1 X 10"^ and 1 X 10"**, respectively (see j|4. ll for discussion). The solid, 
dotted, and dashed curves outline the regions covered by the fa 2 Ms CDF-S, 
GOODS-S, and MUSIC surveys, respectively. The cross near the center in- 
dicates the average aim point of the CDF-S. 

primary counterparts. These counterparts are either selected 
from the SIMPLE catalog (less-accurate positions), VLA cat- 
alog (low background source density), or are relatively bright, 
so that the likelihood ratio could still exceed the threshold at 
large ONIR-X-ray separations (see ^2.2\ . Figure|5]shows the 
distributions of the positional offsets in four bins of differ- 
ent X-ray positional uncertainties. The positional offsets are 
generally smaller than or comparable to the « 85% X-ray po- 
sitional uncertainties (dashed lines in Fig.|5]) for most sources. 

2.4. Comparison with the Error-Circle Matching Method 

To check the efficiency of the likelihood-ratio technique, 
we compared the matching results above to those derived 
from the simple error-circle method described in §3.3.1 of 
LOS. 2^ For the WFI /J-band, LOS gives 347 matches with an 
expected false-match probability of 12% (305 expected true 
matches), while here we have 344 matches with a false-match 
probability of 9% (312 expected true matches). As expected, 
the likelihood-ratio technique provides more reliable matches 
(« 7 more sources using our approach). Between the 347 
counterparts from the error-circle method and the 344 coun- 
terparts from the likelihood-ratio method (if a source has mul- 
tiple counterparts, we chose the one with the highest reliabil- 
ity parameter), 1 1 are different matches, six of which have a 

In LOS, we searched for the closest optical source within a matching 
radius of 1.5^ A^ + cr^^ where Ax is the 85% confidence-level X-ray po- 
sitional eiTor. 
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Positional Offset (orcsec) 

Fig. 5. — Histograms showing the distributions of the positional offset be- 
tween an X-ray source and its primary counterpart. X-ray sources are divided 
into four bins based on their 85% positional uncertainties (Ax) given in 
L08: 0"-0"5, 0"5-l"0, 1"0-1"5, and l"5-2"0. The vertical dashed lines 
indicate the median positional uncertainties (L08) for X-ray sources in each 
bin. The positional offsets are generally smaller than or comparable to the 
X-ray Ri 85% positional uncertainties for most of the sources. 

likelihood-ratio counterpart with magnitude of 7? > 25.5. 

We studied the relation between the false-match probability 
and catalog depth by setting different /?-band magnitude lim- 
its. The results are shown in Figure|6] which indicate that the 
likelihood-ratio method outperforms the error-circle method 
at faint magnitudes, and this is consistent with the expecta- 
tion that the likelihood-ratio method can guard against false 
matches to faint sources. Indeed, the six different counter- 
parts with R> 25.5 derived from the two methods contribute 
to 20% of the total population of 30 > 25.5 count erparts for 
the entire X-ray sample (see also lBrusa et al.ll2b09h . 

From the cumulative distributions of false-match probabil- 
ities shown in Figure |6] we calculated the differential val- 
ues at several /?-band magnitudes. For the likelihood-ratio 
(error-circle) method, the false-match probabilities are « 7%, 
16%, 20%, and 24% (« 7%, 21%, 2S%, and 37%) atR = 24, 
25, 26, and 27, respectively, showing the challenge of deriv- 
ing secure matches at faint optical magnitudes, even with the 
likelihood-ratio technique. 

As a further test, we performed the matching with the 
deeper GOODS-S z-band catalog using the error-circle 
method. There are 266 matches and a 20% false-match prob- 
ability, while the likelihood-ratio method yields 259 matches 
and a 9% false-match probability. 14 matches are different 
from the error-circle results. Ten of these different matches 
have z > 25.5, contributing to « 30% of all the counterparts 
with z > 25.5. Note here that the matching results and the 
false-match probability of the error-circle method depend on 
the choice of the matching radius. Using the same matching 
radius as that for the WFI 7? -band might not be optimal for the 
GOODS-S z-band. It is laborious to determine the optimal 
matching radius for the error-circle method, i.e., the matching 
radius that provides a large number of counterparts without 
introducing too many false matches.^"* On the other hand, 
the likelihood-ratio method can perform a similar optimiza- 
tion (in terms of the completeness and reliability parameters) 

^"^ We tried several other matching radii and none of them can provide bet- 
ter results than the likelihood-ratio method. For example, a matching radius 

of l.l^A^ +0-2 (1.3^A^ + ct2) resulted in 252 (261) matches with a 13% 

(16%) false-match probability. 
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Fig. 6. — Cumulative false-match probability as a function of the 
R-band limiting magnitude for the error-circle method (dotted curve) and 
the likelihood-ratio method (solid curve). The false-match prob abilities are 
estimated using the bootstrapping technique described in j)2.3l EiTor bars 
indicate the errors in the mean values of the false-match probabilities at dif- 
ferent magnitude limits. The corresponding differential false-match proba- 
bilities for the likelihood-ratio (eiTor-circle) method are ft: 7%, 16%, 20%, 
and 24% (si 7%, 21%, 28%, and 37%) at R = 24, 25, 26, and 27, respectively. 
The likelihood-ratio method yields more reliable matches at faint magnitudes. 
The error-circle method appears to have a generally smaller false-match prob- 
ability at bright magnitudes {R < 23); however, it also gives a smaller number 
of matches compared to the hkelihood-ratio method, and the numbers of ex- 
pected true matches from these two methods are comparable. 

automatically and more carefuUy, and is thus more robust. 

We conclude that the likelihood-ratio technique can provide 
much more reliable matches than the error-circle method at 
faint optical magnitudes (e.g., R > 25). The resulting matches 
may differ by « 20% in this regime, where we expect to find 
some of the most interesting objects in deep X-ray surveys, 
such as heavily obscured AGNs and AGNs at very high red- 
shifts. When matching to shallower catalogs with low source 
densities, the likelihood-ratio and error-circle methods will 
produce similar results with competitive reliabilities; how- 
ever, the multiple candidates and reliability parameters pro- 
vided by the likelihood-ratio matching are still useful and con- 
venient when considering multiwavelength identifications. 

3. PHOTOMETRIC REDSfflFTS 

The CDF-S area has been targeted intensively by photo- 
metric and spectroscopic surveys, allowing us to determine 
robustly the photometric redshifts (hereafter, photo-z's) of the 
X-ray sources. Compared t o previous studies of the photo-z's 
of the optical sourc es je^g.. lWolfetani2004.) or X -ray sources 
(e.g., IZheng et"aLl |2004|) in the CDF-S region, our photo-z 
work here has the following features and advantages: (1) We 
are utilizing the X-ray sources in the current « 2 Ms Chandra 
catalog, with the most reliable counterpart information avail- 
able derived using the likelihood-ratio matching method. (2) 
The most recent UV-to-IR deep surveys provide an unprece- 
dented multiwavelength data set for photo-z calculation; the 
Spitzer, GEMS HST, and GALEX UV data have been used for 
the first time in the photo-z calculation for the CDF-S X-ray 
sources. As we are dealing with a relatively small sample of 
sources, we were able to carry out source deblending care- 
fully and obtain the best-ever photometric data for them (see 
^3.11 ). The likelihood-ratio method is used when matching the 
primary counterpart positions (Columns 4 and 5 in Table |2]i 
to the photometric catalogs, and the false-match probabilities 
are small (< 4%). Details regarding the multiwavelength data 



are described below. (3) A collection of up-to-date secure 
spectroscopic redshifts (hereafter spec-z's) allows better cal- 
ibration of the SED templates and a fair assessment of the 
photo-z quality. 

Of the 446 identified X-ray sources, 220 have secure 
spec-z's. We collected secure spec-z's (> 95% confidence 
level, with sev eral secure spe c tral fe at ures) from the follow- 
ing catalogs: iLe Fev re et af ('2004'), 'Szokoly et al.' (2004^ 
Zheng et a l.1 (120041). iM ignoli et al. (2005), Ravikumar et al] 
( 120071) . IVanzella et all (12008). Pope sso et al.i (2009), and Sil- 
verman et al. (2010, in preparation). A matching radius of 
0.5" was used when matching the primary counterparts to the 
redshift catalogs. A bootstrapping test similar to that in §2.3 
indicates a false-match probability of 0.8% for these spec-z's. 
We refer to these 220 sources with secure spec-z's as the spec- 
troscopic sample. Of the other 226 sources, six are identified 
as stars (see details below in ^3.2b . and we refer to the remain- 
ing 220 sources as the photometric sample. Both the spectro- 
scopic and photometric samples have 220 X-ray sources, and 
their photo-z's have been computed. According to the X-ray 
and optical property analysis in Bauer et al. (2004) and Xue et 
al. (2010, in preparation), more than f» 75% of these sources 
are AGNs (see also §3.3), while the others appear to be star- 
burst or normal galaxies. 

3.1. Multiwavelength Photometric Data 

We constructed a photometric catalog for the 446 identified 
X-ray sources including up to 42 bands of UV-to-IR data (up 
to 35 bands are used in the photo-z calculation for a given 
source; see below), derived from a number of publicly avail- 
able catalogs. For a source undetected in a given band, if the 
non-detection is not caused by the possible blending with a 
nearby source, we estimated its 1 a magnitude upper limit, 
using the backgroun d rms maps generated by S EXTRACTOR 
dBertin & A rnouts 1996). The 1 a instead of the more conser- 
vative 3 cr upper limits were calculated because these upper 
limits are used later in the photo-z calculation (see ^3.21 ). Up- 
per limits are not available for the COMBO- 17 and MUSIC 
bands, and we calculated only the z-band upper limits among 
the HST ACS bands. We prioritized the catalogs as an indica- 
tion of our preference for using the data from these catalogs, 
with a highest priority of 5 and lowest 1 . The priorities are 
determined considering the spatial resolution and number of 
bands available (for the consistency in colors, e.g., the MU- 
SIC f/-band has a higher priority than the VIMOS f/-band). 
Details on the catalogs are given below. 

1. The GALEX UV catalog from GALEX Data Release 4, 
which includes the neai--UV (NUV, Aeff = 2267 A) and 
far-UV (FUV, Aeff = 1516 A) data.^*^ The 5 a limiting 
AB magnitudes for the NUV and FUV bands are 26.0 
and 26.4, respectively. This catalog has a priority of 1. 

2. The GOODS-S deep VIMOS f/-band catalog 
(iNonino et al.ll2009l) . with a 1 cr limiting AB magnitude 
of 29.8. This catalog has a priority of 1. 

3. The COMBO- 17 op tical catalog (U B, V, R , /, plus 12 
narrow-band filters: IWolf et aljr2004 llool . The 10 cr 
limiting AB magnitudes for the U, B, V, R, and / bands 

We cannot apply the likelihood-ratio method described in §2.2 to match 
the primary counterparts to the optical sources with secure redshifts, as these 
sources are not from well-defined and flux-limited catalogs. 
^'' Ihttp : / /galex . stsci . eciu/GR4 /| . 
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are 24.7, 25.0, 24.5, 25.6, and 23.0, respectively. This 
catalog has a priority of 3. 

4. The MUSYC fiV/^-dete cted optical catalog (t/, B, o3, 
y, R, I, and z filters; Gawiser et"aLl 12006 '). and the 
MUSYC /T-detected^ptical/NIR catalog {B, V, R, /, 
z, 7, and K filters; .Taylor et all 120091) . The 5 a Hmit- 
ing AB magnitudes for the J and K bands are 23. 1 and 

22.4, respectively. These catalogs have priorities of 2. 

5. The GOODS-S MUSIC optical/NIR/IR catalog (f/ss, 
J/38, f/, F435W, F606 W, F115 W, F850LP, J, H, K, 
and IRAC filters; GrazianetalJ|2006). The 90% Hm- 
iting AB magnitudes for the J, H, and K bands are 

24.5, 24.3, and 23.8, respectively The {/-band filter is 
the same as that used in the GOODS-S deep VIMOS 
t/-band catalog. This catalog has a priority of 3. 

6. The WFI R-hmd catalog (see §2 of iGiavahsco et all 
|2004), with a 5 limiting AB magnitude of 27.3. The 
/?-band filter is the same as that used in the COMBO- 17 
catalog. This catalog has a priority of 2. 

7. The GOODS-S HST ACS B, V, i, and z-band (F435W, 
F606W, F115W, and F850LP filters) optical catalog, 
with 5 (T limiting AB magnitudes of 29.1, 29.1, 28.5 
and 28.2, respectively. This catalog has a priority of 5. 

8. The GEMS V and z-band (F606W and F850LP) cata- 
log dCaldwell et al.l2008l) .^'^ with 5 a hmiting AB mag- 
nitudes of 28.5 and 27.3, respectively. This catalog has 
a priority of 4. 

9. The SIMPLE IRA C 3.6 ^m, 4.5 ^m , 5.8 ^m, and 
8.0 ^im IR catalog (iDamen et al.ll201d) . with 5 a Hm- 
iting AB magnitudes of 23.8, 23.7, 22.0, and 21.8, re- 
spectively. This catalog has a priority of 1 . 

To collect photometric data from these catalogs for the 
X-ray sources, we matched the primary counterparts of the 
X-ray sources to the UV, optical/NIR, and IR catalogs us- 
ing the same likelihood-ratio method described in 32.21 The 
bootstrapping exercises gave small false-match probabilities 
(^S 4%). The aperture-corrected total magnitudes were ex- 
tracted from the catalogs for each matched source. Figure [T] 
shows the SED of source 312 as an example of the resulting 
multiwavelength photometric data. 

Not all the data are used in the photo-z calculation. We 
ignored the data if the errors are larger than one magnitude. 
When detections in the same filter are present in different cat- 
alogs, we chose the data from the catalog with the highest 
priority. When a source has both MUSYC and COMBO-17 
detections, the MUSYC optical data (except for the o3 and z 
bands) were ignored, as there are a larger number of filters 
in the COMBO-17 data. When a source has both MUSYC 
and MUSIC NIR detections, we neglected the MUSYC NIR 
data. Therefore, any individual source can have up to 35 fil- 
ters of photometric data that are used for deriving the photo- 
z. We prefer the IRAC data in the GOODS-S MUSIC cata- 
log because the IRAC measurements are less affected by pos- 
sible source-blending problems. In cases when photometric 
data in two similar filters (e.g., GOODS-S MUSIC iZ-band 
and MUSYC t/-band) differ by more than 0.5 magnitudes 

See |http : / /www .mpia-hd .mpg . de /GEMS /gems ■ htm for the 
V-band data. 



(which may be caused by source variability or blending), we 
used only the data with the higher priority. For about ten 
GEMS sources, the z and V-band magnitudes are > 1 mag- 
nitude fainter than the MUSYC z and V-band data. We de- 
rived aperture photometry (with a 0.7" aperture radius) re- 
sults using the GEMS science images, and they agree with 
the MUSYC photometry data. These sources appear to have 
some low-surface-brightness emission, and the GEMS cata- 
log underestimates their brightness. Therefore we ignored the 
GEMS data for these sources. Upper-limit information was 
only used when there is no detection in a given band and all 
other similar bands. 

As we already know the presence of the X-ray sources and 
their primary counterparts, we can search for their detections 
in the optical/NIR images at low significance levels with less 
worry about false detections. Therefore, we ran SEXTRAC- 
TOR on the WFI R-hanA and MUSYC optical/NIR images 
with various detection significances down to 1 a; we did not 
search the GOODS-S and GEMS optical images because the 
corresponding catalogs are already very complete. On aver- 
age, only a few additional detections were found in each band, 
and we visually examined the images to verify their validity. 

To reduce the effects of source blending, we inspected each 
source and corrected its photometric data if these are affected 
by blending. About 10% of the sources in the IRAC images 
and a few percent in the optical/NIR are blended with adjacent 
sources, and they are usually resolved in the deep HST im- 
ages, i.e., in the GOODS-S or GEMS z band. We deblended 
these sources by calculating the flux ratios using the pixel val- 
ues at the source positions. In some cases when there is only 
mild blending in the IRAC images or the blended sources are 
not detected in the z band, we manually extracted the aperture 
fluxes by centering on the correct source positions (positions 
of the primary counterparts) and using an aperture radius of 
1.5"; the fluxes were then aperture corrected to give the to- 
tal flux. For the one case where the primary counterpart is 
blended (source 22), we centered on the X-ray position to de- 
rive the IRAC fluxes. We ignored the UV data if the source 
suffers from blending, because the UV image resolution is not 
as good as those for the optical/NIR/IR images. 

The final photometric catalog is presented in Table |4]^^ 
Rows 2-463 of Table |4] list the source X-ray IDs (Column 1) 
as well as the photometric data, data errors, and data flags 
in each band: the GALEX bands (Columns 2-7), VIMOS U 
band (Columns 8-10), COMBO-17 bands (Columns 11-61), 
MUSYC bands (Columns 62-103), MUSIC bands (Columns 
104-145), WFI R band (Columns 146-148), GOODS- 
S ACS bands (Columns 149-160), GEMS ACS bands 
(Columns 161-166), SIMPLE IRAC bands (Columns 167- 
178), and VLA 1.4 GHz band (Columns 179-181). The radio 
data are used only in the X-ray source identification process, 
and are not used in the photo-z calculation. All the photo- 
metric data are given in AB magnitudes, and are set to "99" 
if not available; upper limits are indicated by the "-" sign. 
A value of "1" for the data flag indicates that the data point 
was used in the photo-z calculation, "0" means that the data 
point was not used or not available, and "-1" means that the 
data point is probably problematic, either being blended or 
disagreeing with other data. Row 1 of Table |4] provides, for 
the identified extragalactic sources, the median magnitude (in 
the magnitude column), normalized median absolute devia- 

A FITS file containing Table B] (without the fi rst row) is av ailable at 
|http : / /www ■ astro .psu ■ edu/~niel/cdf s/cdfs -Chandra ■ html| 
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Fig. 7. — The observed median SED for the CDF-S X-ray sources and 
its best-fit SED template (a hybrid type 2 AGN template). Data in the same 
filter are combined (e.g., GOODS-S z and GEMS z, MUSIC U and VIMOS 
U). The NUV, FUV, and MUSIC H bands are not shown as the median mag- 
nitudes were calculated using only the detected sources in these bands, and 
thus the SED consists of 39 bands of data. The MUSIC J, K, and GOODS- 
S B and j-band data only represent X-ray sources in the GOODS-S region. 
The en'ors for the median magnitudes are estimated using NMAD / \/N,„^j . 
The limiting magnitudes for some of the deep optical-to-IR bands are also 
shown. This fit should not be interpreted physically as sources across a broad 
range of redshift are included. 

tion (NMAD; in the magnitude-error column), and number 
of X-ray sources used to derive the median magnitude (A^med; 
in the data flag column) for each band.-^ The median mag- 
nitude in a given band was calculated based on the magni- 
tudes of all the 440 sources (after excluding the 16 uniden- 
tified sources and 6 stars), and for undetected sources, the 
limiting magnitude of that band was used.^*' There are some 
exceptions: for the GOODS-S and MUSIC catalogs, we used 
only X-ray sources within the coverage of these catalogs; for 
the NUV, FUV, and VLA 1.4 GHz bands, in which less than 
half of the X-ray sources are detected, and for the MUSIC 
U and H bands, which cover only a portion of the MUSIC 
field, we used only detected sources to calculate the median 
magnitudes. The data in Row 1 indicate how bright a typi- 
cal CDF-S X-ray source is in each band; note that the median 
values in some bands do not sample the entire CDF-S X-ray 
source population. The SED derived from the median mag- 
nitudes and its best-fit template are shown in Figure [T] (see 
details abo ut the SED template library and photo-z fitting pro- 
cedure in § 13.2b . This fit should not be interpreted physically 
as sources across a broad range of redshift are included. 

FigurelSt shows the distributions of the number of filters in 
which the source is detected, and the number of filters used 
in the photo-z calculation (including upper limit data). The 
maximum values of these two filter numbers are both 35, and 
the minimum values are 3 and 7, respectively. We also show 
in Figure [8]5 that brighter /?-band sources tend to have better 
multiwavelength coverage with a larger number of filters.^' 

The normalized median absolute deviation is a robust measure of the 
spread of the magnitudes (cr) in a given band, defined as NMAD = 1.48 X 
median(|m(AB)-median(m(AB))|) (e.g.. lMaronna et all2066y 

The median magnitude and normalized median absolute deviation ai'e 
stable parameters in general, and are not significantly affected by the choice 
of using the limiting magnitude or the individual magnitude upper limits for 
undetected sources, as long as more than half of the sources are detected. 

' For some sources that do not have a WFI R-band detection, we used the 
MUSYC i?-band magnitudes if available. The MUSYC survey searched for 
sources in the B, V, R, and K bands and thus provides additional detections. 



Sources in the spectroscopic sample are usually much brighter 
and have more photometric data than the other sources. 

3.2. Photometric Redshift Fitting 

We use the Zu rich Extragalactic Baye sian Redshift An- 
alyzer (ZEBRA; iFeldmann et al.l 12006 *) to compute the 
photo-z's. ZEBRA provides a maximum-likelihood estimate 
for the photo-z's of individual galaxies or AGNs, combining 
together some useful features such as automatic corrections of 
the SED templates to best represent the SEDs of real sources 
at different redshifts (the training procedure). Briefly, for ev- 
ery input SED template in each user-specified redshift bin, 
ZEBRA generates a new SED template that simultaneously 
minimizes the values for all the training sources (with 
spec-z's available) that use this input template as the best- 
fit template. The procedure is iterated several times, select- 
ing the best-fit template from also the previous trained tem- 
plates from the secon d iteration onward (see more details in 
IFeldmann et alj|2b06l) . The template correction/training pro- 
cess substantially reduces template mismatches and improves 
the photo-z accuracy. In the template-training mode, we used 
the spectroscopic sample as a training set to improve the SED 
templates. 

To construct a suitable library of S ED templates, we fol - 
lowed a procedure similar to that in ISalvato et al.l (|2009|) . 
which calculated photo-z's for the XMM-Newton COSMOS 
sources. First, we ran ZEBRA on the spectroscopic sample 
of 220 sources, fixing the redshifts to their spec-z's and using 
the 259 PEGASE galaxy templates that cover different galaxy 
types (elliptical, spiral, and starburst) and span a wide range of 
star-f ormation history and intrinsic extinction (iGrazian et alJ 
l2006l) . Of these templates, 77 were used by the 220 spec- 
troscopic sources. Second, we repeated the above p rocedure 
using the 10 AGN templates of iPolletta et alj (l2007h instead, 
which include a variety of empirical quasar and Seyfert tem- 
plates. Eight out of the 10 AGN templates were used, includ- 
ing three type 1 and five type 2 AGNs. Third, to account for 
contributions from both the nuclear emission and host-galaxy 
emission in AGNs, we constructed a series of hybrid tem- 
plates with the five most frequently used galaxy templates and 
four most frequently used AGN templates. For each pair of 
AGN and galaxy templates, the templates were normalized by 
the total 3000^000 A fluxes, and nine hybrids with varying 
galaxy/AGN ratios (from 90:10 to 10:90) were created. This 
procedure produced a set of 265 templates that consists of 180 
hybrid templates, 77 galaxy templates, and 8 AGN templates. 
W e chose a miich larg er sample of SED templates than that 
in lSalvato et alj (120091) . which has only 30 templates, so that 
we are less affected by selection effects when estimating the 
photo-z quality using the spectroscopic sample (see ^33] be- 
low). We have tested several smaller SED libraries, which 
provide little improvement in the photo-z quality for the spec- 
troscopic sample, whereas limited SED templates will proba- 
bly hurt the performance of ZEBRA on the photometric sam- 
ple as they might not represent the typical SEDs of the sources 
without spec-z's. 

When using ZEBRA to fit the SEDs, we set most of the 
parameters to their default values. The filter responses and 
galaxy templates were resampled to a uniform wavelength 
mesh with a resolution of 4 A between 400 A and 10 /im. 
To minimize possible spurious boundary effects, we ran the 
template-correction mode twice with two different sets of 
redshift bins and then averaged the resulting improved tem- 
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Fig. 8. — (a) Distributions of the number of filters used in tlie plioto-z calculation (solid histogram) and the number of filters in which a source is detected 
(dashed histogram). The difference between these two histograms indicates the upper-limit information used in the calculation. The maximum number of filters 
used in the photo-z calculation is 35, and the minimum number is 7. The 16 sources without identifications have a filter number of 0. The apparent bimodalities 
of the distributions are due to the inclusion or exclusion of the COMBO- 17 data, which gives a difference of 17 in the number of filters, (h) Number of filters 
used in the photo-z calculation vs. the WFI/MUSYC S-band magnitude for sources with spec-z's (filled dots) and sources without spec-z's (open dots). For 
62 sources that do not have R-band data, their R-band upper limits are shown. Sources in the spectroscopic sample are typically much brighter and have more 
photometric data than the other sources. A few bright (e.g., R < 22) sources have relatively small numbers of filters, because their COMBO-17 data are not 
available (not detected due to source blending). Also note that for very faint sources, e.g., R > 26, most of the data available are upper limits. 



plates to create a final SED library. For the input photometric 
catalog, we employed a minimum magnitude error of 0.06 
for each filter, which effectively avoids cases where a few 
data points with unrealistically small errors have very large 
weights and thus drive the overall fitting results. We ap- 
plied a luminosity prior of Mg < -18, which forces a lower 
limit on the absolute B-band magnitude. There are only two 
sources in the spectroscopic sample with fainter magnitudes. 
Applying such a prior leads to higher photometric redshifts 
for such faint sources, but tests show that it reduces the per- 
centage of catastrophic photo-z failures (outliers in a spec-z 
vs. photo-z plot) considerably. ZEBRA does not take upper 
limits in the input data. To circumvent this problem, we used 
a "workaround" by treating the 1 a upper limits as actual de- 
tections and setting the flux ranges to be from to 2 cr (R. 
Feldmann 2009, private communication). These data points 
generally have a small weight in the SED fitting because of 
the large error bars, but they do help ZEBRA to determine the 
correct location of the Lyman break and thus give a better es- 
timate of the photo-z, especially for sources with upper limits 
only in the NUV/FUV, and sources with detections only in the 
NIR/IR. 

For each source, the ZEBRA SED fitting procedure gen- 
erates a two-dimensional probability distribution in template 
and redshift space; the probabilities indicate the quality of 
fit based on the maximum-likelihood analysis. The best fit 
with the highest probability (corresponding to the lowest 
value) is chosen as the final photo-z and SED template for 
each source. In the photometric sample, « 55% of the sources 
having multiple local maximum probabilities. As a supple- 
ment, we compare these local maximum probabilities, select 
the fit with the second highest probability (corresponding to 
the second lowest value), and derive an alternative photo-z, 
if the second highest probability is within two orders of mag- 
nitude of the highest value (w 40% of the sources with mul- 
tiple local maximum probabilities). If available, alternative 
photo-z's are provided for the photometric sample in Table |5] 
(see below). 

Stars in the sample were identified by visual inspection of 
the GOODS-S and GEMS HST images and also the SED 



fitting results (the AGN/galaxy templates we used gener- 
ally cannot fit a star SED and there is a strong optical ex- 
cess). Every source was visually examined, though we fo- 
cused on sources identified as stars by the COMBO-17 and 
GOODS-S MUSIC classifications, and sources with low fuU- 
band X-ray-to-z-band flux ratios. There are six stars identi- 
fied, with X-ray source IDs 36, 172, 206, 352, 400, and 420; 
all of them have flux ratios (Fo.s-g.o kev/F,) smaller than 0.01. 

3.3. Photometric Redshift Results 

Photo-z's were obtained for the 440 X-ray sources in the 
spectroscopic and photometric samples. They are listed in 
Table |5] with the details of the columns given below. 

• Column (1): the X-ray source ID number 

• Column (2): the spec-z, set to "-1" if not available. 

• Column (3): the photo-z derived using ZEBRA. Set to 
"0" for the six stars and "-1" for the 16 sources without 
identifications. 

• Columns (4) and (5): the 1 a confidence interval (lower 
and upper bounds) of the photo-z. We caution that the 
error estimate is probably not be a good indicator of the 
accuracy of the photo-z (see 33.4l for more details). 

• Columns (6) and (7): the number of filters in which the 
source is detected, and the number of filters used in the 
photo-z calculation. The difference between these two 
numbers is the number of upper hmits used; see Table|4] 
and Figure[8]for more details. 

• Column (8): the reduced x^ value for the best-fit re- 
sult, Xdof' indicating the quality of fit. A value that de- 
viates significantly from unity is generally caused by 
disagreement between the photometric data in different 
catalogs, which is often likely the result of AGN vari- 
ability. 

• Column (9): the alternative photo-z, set to "-1" if not 
available. There are 49 sources having an alternative 
photo-z. 
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• Column (10): the reduced ^ value for the alternative 

photO-Z, Xdof.alf 

• Column (11): a flag indicating whether the X-ray 
source is detected in the optical. A value of "0" means 
no optical detection, in which case the photo-z was 
calculated using only the IR/NIR data and the optical 
upper-limit information (and is probably not very reli- 
able). There are 20 sources in the photometric sample 
that do not have optical detections. 

• Column (12): S pec-z reference numbe r; numbers 1 to 
8 correspond to 'Le Fevre et al.' (2004*), 'Szokolv et al.' 
(2004), Zhenaetal. (2004), MianoU et al. (2005), 
Ravikum ar et al.l (|2007|) . iVanzeUa et all (l2008l) . 
Po pesso et al.1 (l2009h . and Silverman et al. (2010, in 
preparation), respectively. 

The highest spec-z for a CDF-S X-ray source is 3.70. How- 
ever, the derived photo-z's show that ten sources in the pho- 
tometric sample have Zphoto = 4.3-7.6. These sources are gen- 
erally not detected in the optical, or only have a faint optical 
counterpart. They are detected in 4-10 photometric bands, 
and the resulting photo-z's are probably not very reliable. 
In fact, eight of these objects have an alternative photo-z 
of 2.7-4.4 as listed in Table |5] (for comparison, only 41 of 
the remaining 210 sources in the photometric sample have 
an alternative photo-z). In Figure |9^, we show the SEDs 
along with the best and alternative fitting templates of the 
five sources with Zphoto > 5 (source 95 does not have an al- 
ternative photo-z). The SEDs and fitting templates for the 
five Zphoto « 4-5 sources (XlDs 7, 21, 191, 404, and 439) 
are qualitatively similar to those shown in Figure |9^. Given 
the limited SED data available, both the best and alterna- 
tive fitting results appear acceptable. We estimated the in- 
trinsic hydrogen column densities (A^h; see the end of this 
section for details of A^h estimation) for these sources. Two 
sources (XIDs 7 and 439) have estimated A^h values smaller 
than 6 x 10^^ cm~^, and another seven sources have A^h val- 
ues between 1 x lO^-' and 8 x lO^-' cm"^. Source 261 has a 
very hard X-ray spectrum, which requires a Compton-thick 
column density (A^h ~ 1.9 x 10^^ cm"^; an outlier compared 
to the A^H values of all the other X-ray sources) to produce the 
observed X-ray counts at Zphoto = 5.2, and is thus more likely 
located at its alternative redshift of Zait = 2.7, which only re- 
quires A^H ~ 8 X 10^"' cm~^. Given the properties above and 
the scarcity of z > 4 AGNs in the Chandra deep fields, e.g., 
only three such object s are identified spectroscopically in the 
^ 2 Ms CDF-N (e.g.. I^ddington et al.l ll999t IVignali et al.1 
I200I iBareer et al.ll2003h ^^7 we consider that at least half of 
the ten sources with high photo-z's are probably at interme- 
diate redshifts and the alternative photo-z's are preferable for 
them. In the following analysis, we will still use the best-fit 
p hoto-z's for all the sources. 

iRodighiero et al.l (l2007h estimated photo-z's for a sample 
of optically undetected objects, including three of our X- 
ray sources (XIDs 289, 306, and 322, with photo-z's w 4, 
w 4, and w 8). Their estimated photo-z for source 306 is 
more consistent with its alternative photo-z in our catalog 

^2 One X-ray source (XID 618) in the lGiacconi et al.l <2002r ) catalog for the 
1 Ms CDF-S has a spec-z of 4.759 (iVanzella et aLi'lOOff; Fontanot et al.' 

mThe source was detected only by SEXTRACTOR in Giacconi et al. 
, and was not detected in the Alexander et al.i 12003,) catalog for the 
1 Ms CDF-S or the LOS catalog for the Ri 2 Ms CDF-S. 



(see Fig.|9^). We note that lMainieri et al.l ( l2005h also derived 
a single solution for the photo-z of source 95 (see Fig. |9^), 
though with a higher value, Zphoto = 7.9; this source may be 
a z > 5 AGN candidate. There are also at least two z > 5 
AGN candidates among the unidentified X-ray sources (see 
^4.1b . K oeke moer et al., (2004) studied seven extreme X-ray- 
to-optical ratio sources (EXOs) in the « 1 Ms CDF-S. These 
sources are all included in the L08 catalog, with source IDs of 
95, 98, 133, 189, 235, 290, and 404. Sources 95 and 404 have 
photo-z's of 5.7 and 4.7 in our catalog, respectively, while the 
other five sources have photo-z's ranging from 2.2 to 3.2. It 
thus appears that most of the EXOs are at moderate redshifts 
(e.g..l Mainieri et al.ll2005l) . 

The distribution of all the derived photo-z's is shown in 
Figure [TOk. The median photo-z for all the sources is 1.3, 
and it is 2.0 for just the photometric sample. For compari- 
son, the median spec-z is only 0.7 for the spectroscopic sam- 
ple. Between z = 1 and z = 4, there are 179 sources in 
the photometric sample and only 84 sources in the spectro- 
scopic sample. It is reasonable that sources without spectro- 
scopic information are in general at higher redshifts. We cal- 
culated the rest-frame 0.5-8.0 keV luminosities for the X-ray 
sources using the spec-z's and photo-z's. The intrinsic ab- 
sorption was corrected for by assuming an intrinsic X-ray 
power-law photon index of F = 1.8 and using XSPEC (Ver- 
sion 12.5.0; Arnaud 1996) to derive the appropriate intrin- 
sic hydrogen column density (A^h) that produces the observed 
X-ray band ratio (defined as the ratio of count rates between 
the hard and soft bands). A plot of luminosity versus redshift 
is shown in Figure [TOb . Sources in the photometric sample 
have larger luminosities than spec-z sources on average. We 
selected A GNs from the X-ray sources following the crite- 
ria used by iBauer et all (|2004 . Sources with 0.5-8.0 keV 
X-ray luminosities Lo.5-8.0 kev.rest > 3 x 10"^^ erg s"', effective 
photon indices Feff < 1, or X-ray-to-optical (/?-band) flux ra- 
tios F().5-8.{) kev/-P« > 0.1 are considered to be AGNs. In this 
way, we classified 354 AGNs among the 440 identified X-ray 
sources, 150 in the spectroscopic sample and 204 in the pho- 
tometric sample. It thus appears that the majority of the X- 
ray sources in the photometric sample are AGNs at relatively 
high redshifts. In Figure|9j5, we show the median SED of 330 
AGNs with detections in at least six photometric bands and 
having redshifts z < 4, along with the median SEDs in sev- 
eral different X-ray luminosity bins. The median SEDs in the 
optical are dominated by stellar light from the host galaxies, 
consistent with our SED fitting results (see §3.6). 

3.4. Photometric Redshift Accuracy 

To check the quality of the photo-z's, we compared 
the photo-z's to the spec-z's for the 220 sources in the 
spectroscopic sample (see Figures [TTk . [T2h . and [T3h). 
The catastrophic failures (outliers) are defined as having 
|Az|/(l +Zspec) > 0.15, where Az = Zphoto -Zspec- There are 
3 outliers (1.4%) in the spectroscopic sample. The out- 
liers typically still have acceptable SED fitting results, in- 
cluding small photo-z errors, normal reduced x values, and 
no obvious alternative photo-z choices, and thus cannot be 
discriminated without the spec-z information. Specifically, 
there are two relatively bright outliers in Figure [T3b (sources 
243 and 134, with R « 20.0 and R « 22.8). One of them 

This median spectroscopic redshift may be partially related to the two 
apparent redshift spikes at z = 0.67 and 0.73 in the CDF-S fe.g. JGilli et al.l 
l2003tlSlive"rman et al.l200g) . 
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Fig. 9. — (a) SEDs and best-fit templates (solid curves) for the five X-ray sources with photo-z's greater than 5. Also shown are the alternative fitting templates 
(dashed curves) for four of the sources. These templates are either galaxy templates or hybrid AGN templates. None of the five sources is detected blueward of 
K band (the X'-band magnitude is indicated in each panel if available). The 1 a upper limits on the z, J, and K bands are shown if there are no detections. The 
1 (J upper hmits do not set tight constraints on the actual fluxes, so they may appear below the the SED templates in some cases, (h) The median SEDs of z = 0-A 
AGNs in the CDF-S, derived from the photometric data and best-fit SED templates. The 24 /^m data from the Spitzer Far Infrared Deep Extragalactic Legacy 
Survey (FIDEL; Dickinson et al. 2010, in preparation) were also used. The cyan long-dashed, magenta dotted, black solid, red dash-dotted, and blue dashed 
curves represent median SEDs in different X-ray luminosity bins (as labeled, alo ng with the number of AGNs contributing to each median SED). The grey solid 
curve represents the low-luminosity quasar mean SED from Richar ds et alj 12006 ) with the host-galaxy contribution subfi'acted. The SEDs of CDF-S AGNs are 
generally dominated by host-galaxy light and show heavy extinction in the UV. [See the electronic edition of the Journal for a color version of this figure.^ 




Fig. 10. — (a) Histogram showing the disfiibution of photo-z's. The shaded region shows the distribution of spec-z's. Sources with only photo-z's are 
generally at higher redshifts than spec-z sources, (b) Redshift-luminosity distribution for spec-z sources (filled dots) and photo-z sources (open dots). The solid 
curve represents the limiting 0.5-8.0 keV luminosity at the center of the CDF-S region (LOS), assuming F = 1.8 and Wh = 10"^ cm"^. The dotted line indicates 
io. 5-8.0 kcV.iest = 3 X 10"*^ erg s"' . Some of the photo-z's for the ten objects with Zphoto > 4 are probably not reliable. 



(source 243) is caused by the bright GALEX FUV detection 
(mpuv = 21.1), which should have been suppressed due to the 
Lyman break for a source at z^pec = 1-03; the FUV flux may 
be contributed by a few surrounding faint sources. The other 
one is produced by our luminosity prior of Mb < -18, which 
pushes the photo-z to a larger incorrect value (z,spec = 0.12 and 
Zphoto=0.64). 

We also use two other indicators to evaluate the photo-z 
quality, the normalized median absolute deviation of Az 
(cnmad) and the average absolute scatter (AAS), which are 
defined as 



1 A- /|Az-median(Az) 

ctnmad = 1 .48 X median - 



.spec 



and 



AAS = mean 



|Az| 

1+Zsne 



(6) 



(7) 



cnmad is a robust indicator of the photo-z accuracy af- 
ter the exclusion o f outliers (e.g., ' Maronna et al.l 120061: 
iBrammer et al.ll2008 l). while AAS shows the mean scatter in- 
cluding the contribution from outliers. For our spectroscopic 
sample, (Tnmad = 0.010 and AAS = 0.013. These parameters, 
along with FigurefTTh. indicate that the photo-z's for the spec- 
troscopic sample are very reliable in terms of both the outlier 
percentage and accuracy. We show in Figure[T]the SED fitting 
result for source 312 as an example; the derived photo-z for 
this source is accurate to within 1% [| Az|/(1 +Zspec) = 0.005]. 

However, the above evaluation does not represent the qual- 
ity of the photo-z's for sources in the photometric sample, 
because the SED templates were corrected using the spec-z 
information, and we are thus biased to get optimal fitting re- 
sults for sources with spec-z's. Therefore, we performed a 
"blind test" to derive a more realistic estimate of the over- 
all photo-z accuracy. We randomly picked 3/4 of the sources 
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in the spectroscopic sample, and only the spec-z's of these 
sources were used to train the SED templates. The photo-z's 
of the other 1/4 of the sources are calculated based on the 
improved templates, and these are the blind-test sources as 
there was no prior knowledge utilized about their spec-z's. 
The above process was repeated 20 times to create a statis- 
tically significant sample including 1110 blind-test sources 
(the blind-test sample). Applying the same quality indicators 
above to the blind-test sample, we get an outlier percentage 
of 8.6%, (Tnmad = 0.059, and AAS = 0.065. The results are 
shown in Figures [TTb . [TZb . and \T3h . The photo-z accuracy 
does not appear to degrade at high redshifts (z « 2-4), and the 
accuracy declines at faint /?-band magnitudes (the outlier frac- 
tions for the R > 22.5 and R < 22.5 samples are 11.4% and 
6. 1 %, respectively). Note that among the 1110 sources in the 
blind-test sample, some have the same spec-z, photo-z, and 
R magnitude (same source in the spectroscopic sample, and 
different training samples do not alter the best-fit photo-z), 
and thus cannot be distinguished in FiguresfTTb and[T3b (i.e., 
they appear as the same data point). The effective size of the 
blind-test sample should be around 220, as we only have this 
number of sources in the spectroscopic sample. 

In Figures [T4k and [T4b . we show the distributions of the 
1 (7 photo-z errors estimated from the ZEBRA SED fitting for 
the 440 X-ray sources (see Table |5]l and the blind-test sam- 
ple, respectively. The spectroscopic sources generally have 
smaller photo-z errors than the photometric sources, mainly 
because of their better multiwavelength coverage; the ZE- 
BRA photo-z errors also appear to decrease toward brighter 
magnitudes, for the same reason. Sources with an alternative 
photo-z (blue squares in Fig.[l4h) tend to have a large photo-z 
error, and the outliers (red stars in Figs. [T4k and [l4b ) can- 
not be distinguished based on their photo-z errors. The ZE- 
BRA photo-z errors were calculated in the two-dimensional 
redshift-template space. However, the template plays a minor 
role in the error calculation, unless there are multiple compet- 
itive SED fits with very close values. Therefore, the er- 
rors for most sources were determined using essentially only 
the best-fit template, and are generally underestimates of the 
real errors. We compared the ZEBRA 1 a photo-z errors to 
the real errors (|Az|) using the spectroscopic and blind-test 
samples, and the results are shown in Figure [T4b . The ZE- 
BRA photo-z errors only account for « 30% and « 17% of 
the real errors on average for the spectroscopic and blind-test 
samples, respectively. Similar behavior has been discovered 
for the photo-z errors de rived using other photo-z codes (e.g., 
iHildebrandt eTallllOOSl) . in the sense that the photo-z accu- 
racy is not tightly correlated with the error estimate. A more 
realistic estimate of the photo-z errors for our sources is given 
by the ctnmad parameter, about 1% for the spectroscopic sam- 
ple and 6% for the blind-test sample. 

Although we give an unbiased estimate of the photo-z qual- 
ity using the blind-test sample, we note that there are three 
factors that could affect the results: (1) In the blind-test sam- 
ple, only 3/4 of the spec-z sources are used to train the tem- 
plates, while in the final results, all the spec-z sources are used 
in the training. Thus the photo-z quality should be some- 
what better due to the larger training sample. (2) The qual- 
ity estimate is solely based on the spectroscopic sample, in 
which sources are generally brighter and have better photo- 
metric sampling than the sources in the photometric sample 
(see Fig.[8]3). The photo-z quality for the latter might not be 
as goo d as what we derived here. Indeed, ,Hildebrandt et alj 
( 120081) have shown that the outlier fractions for the latter sam- 



ples are typically larger by at least a factor of two. We per- 
formed another blind test, choosing the 182 spec-z sources 
with WFl /? < 24 as the training sample, and then calculating 
the photo-z's of the other 38 spec-z sources with R ss 24-26. 
The resulting outlier fraction is 10.5% and (Jnmad = 0.067,-^"^ 
which are not greatly increased relative to the 8.6% outlier 
fraction and unmad = 0.059 for the blind-test sample above. 
Note that these 38 faint spec-z sources are still significantly 
brighter than the sources in the photometric sample on av- 
erage (see, e.g.. Figure [Sj)); therefore an outlier fraction of 
« 15-25% is probable for the faint sources (R > 26) in the 
photometric sample. (3) The spec-z's were used when we se- 
lected the SED templates. This should not affect the blind test 
significantly, as we chose a large and relatively complete SED 
library. 

In summary, among the 220 photo-z's in the photometric 
sample, we expect there are « 9% outliers for the « 1 10 rel- 
atively brighter sources (R < 26), and the outlier fraction will 
rise to « 15-25% for the fainter sources (R > 26). The typi- 
cal photo-z accuracy is w 6-7% (ctnmad)- The outlier fraction 
and photo-z accuracy do not appear to have a redshift depen- 
dence (for z ~ 0^). Among the faint R >26 sources, 20 
objects do not have any optical counterparts (see Table |5] in- 
cluding the five Zphoto > 5 sources shown in Fig. |9^), and the 
photo-z's were calculated using only the IR/NIR photometric 
data and the optical upper-limit information. Caution should 
be taken when using these photo-z's. Moreover, the alterna- 
tive photo-z's provided for the 49 sources may occasionally 
be a better choice than the best-fit photo-z's (e.g., for some of 
the Zphoto > 4 sources). 

Recently, the GOODS/VIMOS spectroscopic campaign has 
provided two additio nal secure spec-z's for our X-ray sources 
i lBalestra et ai]|20Tol) . z^pec = 0.740 for source 175 and z^nec = 
2.586 for source 372. Their photo-z's in our catalog are 0.56 
and 2.71, with accuracies of | Az|/(1 +Zspec) ~ 10% and 3.6%, 
respectively. 

3.5. Comparison with Other AGN Photometric-Redshift 
Accuracies 

The photo-z quality of the spectroscopic sample (1.4% 
outlier and ctnmad = 0.01) is comparable to thos e of recent 
photo -z studies fo r the COSMOS op tical (e.g., Ilbertet al] 
120091) and X-ray (ISalvato et alJ l2009h sources. Compared 
to the blind-t est result (46 sources) for the COSMOS 
X-ray sources ( ISalvato et alj |2b09) with 10.9% outliers and 
cnmad ~ 0.023, our result has a slighter smaller fraction 
(8.6%) of outliers but a larger o-nmad (0.059). Note that the 
CDF-S and COSMOS X-ray sources are two different popu- 
lations. The « 2 Ms CDF-S sources have an average Chandra 
effective exposure of ss 1 Ms, while the COSMOS sources 
have an average effective XMM-Newton exposure of « 50 ks. 
In our spectroscopic sample, « 38% of the sources have HST 
i > 22.5, and in the COSMOS spectroscopic sample, only 
« 10% of the sources have Ca nada-France-Hawaii Telescope 
/* > 22.5 (ISalvato et alj|2009h . Our photo-z quality is at least 
comparable to that for the COSMOS X-ray sources, consider- 
ing that we are dealing with a fainter and thus more challeng- 
ing sample. One advantage in the photo-z work for the COS- 
MOS X-ray sources is that a correction was made for optical 

The p hoto -z quality does not appear to degrade from R > 22.5 to i? > 24 
(see Figure fnb ). maybe due to the small sample size of R > 24 sources. 
Therefore the outlier fraction here is comparable to that (11.4%) for the R > 
22.5 sources in the blind-test sample above. 
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Fig. 1 1 . — Comparison of the photo-z's a!fiJ<Spec-z's for (a) the spectroscopic sample and (i>) the Wind-test sample. Thtfapper panels show the distribution 
of Zspcc vs. Zphoto, and the lower panels show Zspcc vs. Az/(1 +Zspec). The sohd lines indicate Zphoto = Zspec, and the dotted lines represent relations of Zphoto = 
Zspec ± 0.15(1 +Zspec), which define the outlier limits. Note that one data point can represent several sources in the blind-test sample (fc). The photo-z's and 
spec-z's agree very well for the spectroscopic sample, while the photo-z's and spec-z's differences are considerably larger for the blind-test sample. The photo-z 
accuracy does not appear to degrade at high redshifts (z « 
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Fig. 12. — The distribution of the photo-z accuracy, Az/(1+Zspec), for (a) the spectroscopic sample and (6) the blind-test sample. 



0.4 



0.6 



0.6 
0.4 

0.2 

n" 

+ -0.0 
\ 

^ -0.2 
-0.4 
-0.6 



If 



20 22 24 

WFI R-band Macnitudc 



26 



0.6 
0.4 



-0.2 
-0.4 



(b) 



If 



20 22 24 

WFI R-band Magnitude 



26 
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14. — The distributions of the photo-z errors derived from the ZEBRA SED fitting for (a) the spectroscopic (filled dots) and photometric (open dots) 
samples and (h) the blind-test sample. The ZEBRA photo-z error is calculated 

Zerr — (flipper — Slower The red Stars represent outliers in the spectroscopic and 
blind-test samples. The blue squares represent sources with an alternative photo-z. In (c), we show the distribution of the ZEBRA photo-z errors in units of the 
real photo-z en'ors (| Az|) for the spectroscopic sample and the blind-test sample. The ZEBRA photo-z eiTors generally underestimate the real errors by a factor 
of ft! 3-6. [See the electronic edition of the Journal for a color version of this figure.] 



17 



AGN va riability, which help s to reduce the outlier fraction and 
cnmad (ISalvato et alJl2009h . Given the nature of the photo- 
metric data in the CDF-S, i.e., multiwavelength observations 
are not separated in well-defined epochs, source variability is 
difficult to take into account. However, we also note that for 
the moderate-luminosity AGNs in our sample, the SEDs are 
often dominated by host-galaxy light (see Fig. |9j3), and thus 
AGN variability may play a smaller role for our sources than 
for the brighter COSMOS AGNs. 

Compared to the photo -z quality for t he 342 X-ray sources 
in the « 1 Ms CDF-S (Zhen g et alj2004l) . with « 15% outliers 
and ctnmad ~ 0.08-0.1,^^ or the ph oto-z quality of th e 429 
X-ray sources in the « 2 Ms CDF-N (iBarger et al.l2003h . with 
94% of the photo-z's accurate to within w 25%, our photo-z 
accuracy estimated based upon the blind-test sample is signif- 
icantly (> 50%) improved. 

3.6. AGN Classification and Best-Fit SED Templates 

The intrinsic hydrogen column density A^h estimated in ^3.31 
(assuming an underlying F = 1.8) can provide a basic AGN 
classification from the X-ray point of view, i.e., "absorbed" 
with A^H > 10-^ cm"^ or "unabsorbed" with A^h < 10^^ cm"^. 
This method has limited diagnostic utility because of the com- 
plexity of AGN X-ray spectra. However, it suggests that 
« 256 (w 72%) of the 354 AGNs in the spectroscopic and 
photometric samples are absorbed, which is in general agree- 
ment with the Rs 60-75% absorbed AGN fraction in the « 1 
Ms CDF- S and 2 Ms CDF-N based on X-ra y spectral anal- 
ysis (e.g..' Bauer eTani2004t llbzzi et al.ll2006l) . 

We investigated the best-fit SED templates selected for the 
440 sources; about 98% of them are improved templates from 
the template-training procedure, indicating the importance of 
the training step. Of the 265 input templates (including their 
trained varieties), 105 are selected as the best-fit SED tem- 
plate for multiple sources and 107 are never used; the most 
frequently used template (a starburst-galaxy template) has 
been chosen for 22 sources. The best-fit SED templates can 
also provide a basic SED-based classification of the X-ray 
sources. Of the 440 source SEDs, 212 are best fit by galaxy 
templates, and the others are fit by AGN or hybrid templates, 
including 84 type 1 AGNs and 144 type 2 AGNs. For sources 
that are detected in only a few or none of the optical bands, 
this classification may be highly uncertain. 

AGNs identified in the X-ray do not completely agree with 
the SED fitting results, as expected based upon disagreements 
between X-ray based and optical-spectroscopy based classi- 
fication schemes (e.g.. Matt 2002; Trouille et al. 2009). For 
example, of the 303 AGNs with detections in at least ten 
UV-to-lR bands, 125 (41.3%) ai-e best fit by galaxy tem- 
plates. Despite the uncertainties in the SED fitting and the 
calculation of X-ray luminosities, some of these sources are 
likel y the X-ray brig ht optically normal galaxies (XBONGs; 
e.g.. IComastri et an i2002). which have high X-ray luminosi- 
ties but do not appear as AGNs in the optical. The na- 
ture of XBONGs remains somewhat mysterious: the prop- 
erties of XBONGs are usually explained by heavy obscura- 
tion covering a large solid angle (e.g.. Matt 2002; Rigby et aL 
2006; Civano et al. 2007 ), dilution by host galaxies (e.g. , 
Moran et alji2002i; iSevergnini et al.ll2003l; iTrump et al.ll2009l; 

We also checked the photo-z quaUty of thel Zheng et alj 12004*) sources 
using 31 spec-z sources that did not have spec-z information in Zheng et al. 
t2004ft . The resulting outlier fraction is 26%, and the outlier-excluded 
accuracy is ctnmad ~ 0.08. 



see also Fig . [9b), or a radiatively inefficient accretio n flow 
(RIAF; e.g.. lYuan & NaravanlllOOi: iTrump et alJl2009h . 

An apparent inconsistency also exists between the X-ray 
(unabsorbed versus absorbed) and SED-based (type 1 ver- 
sus non-type 1, including type 2 and galaxies) classification 
schemes. Figure [15] displays the X-ray band ratio versus red- 
shift for the 354 AGNs. The SED-based classification is color 
coded, while the X-ray classification is indicated by the the- 
oretical band-ratio-redshift track. We study the classifica- 
tions in two X-ray luminosity bins, divided by our median 
X-ray luminosity, Lo.s-s.o kev.rest ~ 2.6 x lO"*"* erg s"'. We 
consider only the 303 AGNs with detections in at least ten 
UV-to-lR bands. In the low-luminosity (high-luminosity) 
bin, 42% (55%) of the type 1 AGNs are unabsorbed, 64% 
(86%) of the type 2 AGNs are absorbed, and 71% (87%) of 
the AGNs fit by galaxy templates are absorbed. For rela- 
tively luminous X-ray sources, the two classification schemes 
tend to agree better, probably due to less fractional contam- 
ination from host-galaxy light. It appears that most of the 
non-type 1 sources (ss 86% for luminous X-ray sources) are 
absorbed, while only about half of the type 1 sources are un- 
absorbed. The comparison between t he X-ray and optical- 
spectroscopy classification schemes in iTrouille et alJ (l2009h 
shows that ss 80% of the broad-line AGNs are unabsorbed and 
« 67% of the non-broad-line AGNs are absorbed. Therefore, 
the X-ray versus optical-spectroscopy classification schemes 
agree better for unabsorbed/type 1 AGNs, compared to the X- 
ray versus SED classification schemes. Note that the uncer- 
tainties in the simple X-ray classification scheme above are 
plausibly responsible for some of the classification discrep- 
ancy; detailed X-ray spectral analysis will provide a better 
X-ray classification and give more insight into the problem of 
X-ray versus SED classification inconsistency. 

4. DISCUSSION 

4.1. Nature of the Unidentified X-ray Sources 

There are 16 unidentified sources, eight of which are prob- 
ably related to a star (one source), galaxy groups and clus- 
ters (four sources), off-nuclear sources (two sources), or any 
of a number of galaxies at the edge of the CDF-S field (one 
source; see Col. 13 of Table [2]). For the other eight X-ray 
sources (flagged with "NolD" in Col. 13 of Table [2 see Fig- 
ure [4] for their locations), six were detected only at a false- 
positive threshold of 1 x 10"^, and two (sources 49 and 345) 
at 1 X 10"^. At a threshold of 1 x 10"^, we expect less than 
0.2 false detections (see §3.2 of L08). We visually checked 
the X-ray images of sources 49 and 345 and confirmed their 
apparent validity. We also investigated the Chandra events 
for these two sources and excluded the possibility of con- 
tamination from short-lived cosmic-ray afterglows. There- 
fore sources 49 and 345 are likely genuine X-ray sources. 
Source 49 is outside of the GOODS-S region. It was de- 
tected in all three X-ray bands (« 117 full-band counts) and 
was also present in th e w 1 Ms CDF-S ( Giacconi et alJ[2002l : 
[Alexander et'ai]|2003h . It is a hard X-ray source with a band 
ratio of 1 .5 (corresponding to an intrinsic column density of 
A^H ~ 6 X 10^^ cm"2 assuming z = 5 and F = 1 .8). Its 0.7-6 keV 
X-ray spectrum shows no obvious emission-line features. The 
X-ray-to-optical (z-band) flux ratio is ^b.s-s.o kev/-Pz ^ 50. 

IS_alvato_et alj (2009) also made a similar comparison of the SED versus 
X-ray classifications for the COSMOS X-ray sources, and found a generally 
better agreement (si 80%), probably due to the fact that COSMOS samples 
brighter X-ray fluxes than the CDF-S. 
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Fig. 15. — X-ray band ratio vs. redshift for the (a) fainter and (b) brighter halves of AGNs, separated by the median X-ray luminosity, Lq. 5-8.0 keV.rest = 
2.6 X 10*^^ erg s"'. Band ratios are defined as the ratio of count rates between the hard and soft bands. The green crosses, blue triangles, and red squares are 
sources whose SEDs are best-fit by type 1 , type 2, and galaxy templates, respectively. Upper limits (lower hmits) on the band ratios are shown when the sources 
are not detected in the hard band (soft band). The en'ors on the band ratios are not shown for clarity; the median error is ft: 0.16. The sohd (dashed) curve and 
the dark (light) shaded region show the expected band ratios for AGNs with power-law photon indices F = 1.8 it 0.5 and intrinsic absorption Na = 10^^ cm"^ 
(Nn = 10-^ cm"^). The X-ray (unabsorbed versus absorbed) and SED-based (type 1 versus non-type 1) classification schemes agree better for luminous AGNs; 
~ 86% of the non-type 1 sources are absorbed, and about half of the type 1 sources are unabsorbed. [See the electronic edition of the Journal for a color version 
of this figure.} 



There appears to be marginal 3.6, 4.5, and 8.0 ftm emission 
from the X-ray position of source 49. However, there is a 
relatively bright [m(3.6 ^m) = 21.2] IR source nearby (2.3" 
away); even if the emission is from source 49, it is completely 
blended with the bright source. There is also a faint VIMOS 
f/-band source {U = 28 .70 ± 0.22) 1.3" away from source 49, 
which is not detected in any of the other bands. It was not con- 
sidered a counterpart by the likelihood-ratio matching. Source 
345 is inside the GOODS-S region and is only detected in the 
soft band with w 6 counts. The detection is likely significant 
because there are five counts landing on a single pixel. Its 
X-ray-to-optical flux ratio is Fo.5-8.0 kev/^i ^ 10. The near- 
est ONIR source to source 345 is w 2.0" away. 

The properti es of sources 49 and 345 are comparable to 
the EXOs in Ko ekemoer et alJ (120041) . except that they do not 
have any clear detections in the IR or NIR bands (which EXOs 
generally have) and source 345 is relatively X-ray weak. Most 
of the EXOs are suggested to be dusty galaxies at moder- 
ate redshifts (z « 2-3), and a small fraction could be AGNs 
at very high redshif ts (z > 5; e.g., iKoekemoer et al. 2004; 
iMainieri et alJ l2005h . Given the apparent non-detection of 
these two sources in the deep IR and NIR images, they are 
unlikely to be dusty galaxies at moderate redshifts (i.e., their 
host galaxies should be detected). Therefore, sources 49 and 
345 are more likely to be AGNs at very high redshifts that 
appear blank in all the current optical-to-lR catalogs due to 
finite-sensitivity effects. Using the luminosities and spec- 
z's or photo-z's of our X-ray sources, and the depths of the 
optical-to-IR catalogs, we estimate that w 45 (« 10%) of 
the sources would appear blank in all these bands if placed 
at z = 5, and « 10 of them would still be detected in the X- 
ray if at z = 5. There is also the possibility that source 345 is 
an isolated neutron star, such as th e so-called "Dim Isolated 
Neutron Stars" (e.g.. lKaplanll2008h . which emit thermal radi- 
ation with temperatures of « 0.1 keV. This is not likely the 
case for source 49 because of its hard X-ray emission. Iso- 
lated neutron stars with faint X-ray fluxes such as source 345 
would appear blank at optical-to-lR wavelengths, because 
of the high (w lO^'-lO'*) X-ray-to-optical flux ratios of these 



sources. Given the X-ray sensitivity of the CDF-S (L0 8) and 
the number counts of isolated neutron stars (e.g., P opov et alj 
2000, 2010), we expect « 0.2 detectable isolated neutron stars 
in the CDF-S field. Very deep NIR imaging or spectroscopic 
observations will help elucidate the nature of sources 49 and 
345. 

The six sources (XIDs 93, 106, 210, 280, 344, and 353) 
detected at a low-significance level of 1 x 10"^ and without 
identifications could be spurious X-ray detections. Visual 
inspection indicates that most of these sources appear only 
marginally above the background. At a detection threshold of 
1 X 10"*, we conservatively estimated that there are « 18 false 
detections (see § 3.2 of LOS), and the t rue number is likely 
« 2-3 times less ([Alexander et alJl2003h . Spurious X-ray de- 
tections could also exist in the identified source sample be- 
cause of chance correlations. We estimate that the average 
possibility of matching a randomly placed 1 x 10~* X-ray 
source to an ONIR source is « 11% , using the same boot- 
strapping technique described in 32. 3P ^ If the six unidentified 
X-ray sources are indeed spurious, we expect less than one 
[6 X 1 1%/(1 - 11%)] spurious X-ray source in our identified 
source sample. Therefore, there are up to « 7 spurious X-ray 
detections, six unidentified and one falsely identified. This 
number is roughly consistent with the expected value (« 2-3 
times less than 18). If any of these six low-significance and 
unidentified sources is a real X-ray detection, it could have a 
similar nature as sources 49 and 345. 

4.2. Future Prospects for Improved Source Identification 

Besides the unidentified X-ray sources, uncertainties are 
still present for the 442 identified sources (excluding the four 
additional likely cases), as there is a w 6.2% false-match 
probability and there are 21 primary counterparts from the 
MUSYC or SIMPLE catalogs with relatively less-accurate 

Essentially, this is an estimate of the false-match probability for X-ray 
sources detected at only the 1 X 10"* significance level. The false-match 
probability is higher than that for the entire X-ray sample (Ri 6.2%), mainly 
because the positional uncertainties of these X-ray sources are relatively 
larger 
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(0.2-0.3") positions. In order to characterize the nature of 
all the X-ray sources in the w 2 Ms CDF-S, efforts should fo- 
cus on obtaining (a) better positions for the counterparts and 
(b) a lower overall false-match probability. Better counter- 
part positions generally require deeper optical or NIR obser- 
vations. Comparing the current deepest optical (GOODS-S 
z-band) and NIR (MUSIC /T-band) catalogs, we consider the 
NIR band to be probably the most beneficial to push deeper, 
because of its higher counterpart recovery rate, X-0 fraction, 
and lower false-match probability (see Table [T]! at the current 
depth. Deep HST WFC3/IR surveys or the forthcoming James 
Webb Space Telescope (JWST) will provide unprecedented 
NIR/IR data for source identification. Note that a deeper opti- 
cal or NIR catalog will generally result in a larger false-match 
probability. 

There are at least two possible approaches toward a smaller 
overall false-match probability: deeper X-ray observations 
and deeper radio observations. The probability distribution of 
angular separations, /(r), is largely controlled by the X-ray 
positional errors (see Eq. 2) except when matching to the 
IR sources. Additional X-ray observations will improve the 
X-ray positional accuracy and thus reduce the false-match 
probability. For example, as the « 2 Ms CDF-S is sched- 
uled to be expanded to w 4 Ms, the factor of « 2 increase 
in counts will reduce the areas of source positional error re- 
gions by « 30% on average, reducing the number of spurious 
matches (at the current flux limit) by about this same factor 
The VLA radio catalog has the highest X-O fraction (« 28%) 
and lowest false-match probability (see Table [TJ. However, 
it also has the lowest counterpart recovery rate (« 20%). To 
examine if the high X-O fraction of the radio catalog is intrin- 
sic or rather just related to the low recovery rate (i.e., limited 
catalog depth), we compared it to the X-O fractions of op- 
tical and NIR catalogs at similar effective depth. We limit 
the GOODS-S z-band catalog and MUSIC K-hand catalog to 
brighter magnitudes (z < 21.0 and K < 19.8) until the recov- 
ery rates of these two catalogs are the same as that of the radio 
catalog (sa 20%). The resulting X-O fractions are only « 13% 
and « 17% for these two catalogs, respectively. This suggests 
the X-ray sources are more strongly associated with the radio 
sources than optical or NIR sources at the current depth of the 
radio catalog, and thus deeper radio observations might well 
provide more reliable matches than the GOODS-S or MUSIC 
catalogs. The current radio catalog has a 5 cr limiting flux of 
w 40 mJV- G iven the VLA 1 .4 GHz number counts estimated 
in lOwen Morrison! (l2008h . there will be w 20 times more 
radio sources at a flux limit of 5 /iJy. This radio source density 
(about 30 sources arcmin"^) is still « 2^ times smaller than 
those of the optical and NIR catalogs, and we expect that such 
deeper radio catalogs can provide more reliable matches and 
thus reduce the overall false-match probability. The flux limit 
of 5 /iJy can be easily achieved by facilities such as the Ex- 
panded Very Large Array (EVLA). The Atacama Large Mil- 
limeter Array (ALMA) at shorter wavelengths may also help 
with the goal of obtaining more reliable counterparts in a sim- 
ilar way. 

5. SUMMARY OF RESULTS 

1. We presented multiwavelength identifications of the 
462 X-ray sources in the main Chandra source cata- 
log for the « 2 Ms CDF-S. The optical-to-radio cata- 
logs used for source matching include the WFI /?-band, 
GOODS-S z-band, GEMS z-band, GOODS-S MUSIC 
K-hmd, MUSYC /T-band, SIMPLE 3.6 /im, and VLA 



1.4 GHz catalogs. The matching results are combined 
to create lists of primary and secondary counterparts. 
We identified reliable counterparts for 442 (95.7%) of 
the X-ray sources, with an expected false-match prob- 
ability of « 6.2%; these are the best identifications for 
sources at these X-ray fluxes. We also associated four 
additional likely counterparts given their small posi- 
tional offsets to the relevant X-ray sources. 

2. The likelihood-ratio method was used for source 
matching, which significantly reduces the false-match 
probability at faint magnitudes compared to the 
error-circle matching method. For example, the 
false-match probabilities at = 26 are « 20% and 
w 28% for the likelihood-ratio and error-circle meth- 
ods, respectively. The resulting matches may differ by 
w 20% for sources with R > 25.5, where we expect to 
find some of the most interesting objects in the deep 
X-ray surveys. The likelihood-ratio matching tech- 
nique can also generate multiple counterpart candidates 
and calculate their relative reliability parameters, which 
are useful when considering multiwavelength identifi- 
cations. 

3. Among the optical-to-radio multiwavelength catalogs 
used for identification, the false-match probabilities 
range from w 1-9%, and the counterpart recovery rates 
range from « 20-80%. Given the current depths of the 
catalogs, the NIR catalogs (e.g., the MUSIC /T-band 
catalog) can provide a large fraction of the counterparts 
for the X-ray sources with a relatively low false-match 
probability. Meanwhile, the radio catalog can provide 
the most secure counterparts with almost no spurious 
matches. Even after limiting the optical/NIR catalogs 
to similar effective depth as the radio catalog, the ra- 
dio sources still have a significantly stronger associa- 
tion with the X-ray sources (measured by the X-O frac- 
tion parameter). Therefore, obtaining deeper radio ob- 
servations will considerably benefit the matching reli- 
ability. Additional X-ray observations will reduce the 
X-ray positional errors and thus also improve the iden- 
tification reliability. 

4. Of the 16 unidentified X-ray sources, eight are probably 
related to a star (one source), galaxy groups and clusters 
(four sources), off-nuclear sources (two sources), or any 
of a number of galaxies at the edge of the CDF-S field 
(one source). Another six sources detected at a WAVDE- 
TECT false-positive threshold of 1 x 10~^ may be spu- 
rious X-ray detections. The remaining two objects are 
likely genuine X-ray sources. Their X-ray-to-optical 
properties are comparable to the EXOs, though without 
any IR or NIR detections. They are likely to be AGNs 
at very high redshifts (z > 5), or one of them (source 
345) could be an isolated neutron star. 

5. We constructed a photometric catalog for the 446 
identified X-ray sources including up to 42 bands of 
UV-to-IR data. We performed source deblending care- 
fully for « 10% of the sources in the IR bands and 
a few percent in the optical and NIR bands, obtain- 
ing the best-ever photometric data for these sources. 
High-quality photo-z's were then derived using the ZE- 
BRA SED fitting code and a set of 265 galaxy, AGN, 
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and galaxy/ AGN hybrid templates. The ZEBRA train- 
ing procedure corrects the SED templates to best rep- 
resent the SEDs of real sources at different redshifts, 
which reduces template mismatches and improves the 
photo-z accuracy. Catalogs of photometric data and 
photo-z's are provided. We collected secure spec-z's 
for 220 of the X-ray sources. It appears that the major- 
ity of the X-ray sources without spec-z's are AGNs at 
relatively high redshifts (z « 1-4). 

6. The photo-z's are extremely accurate for sources with 
spec-z's. There are only 3 outliers out of 220 sources, 
and the photo-z's are accurate to within « 1% (unmad)- 
We performed blind tests to derive a more realistic es- 
timate of the overall photo-z accuracy for sources with- 
out spec-z's. We expect there are w 9% outliers for the 
w 1 10 relatively brighter sources (R < 26), and the out- 
lier fraction will increase to « 15-25% for the « 110 
fainter sources {R > 26). The typical photo-z accuracy 



is w 6-7% ((Tnmad)- The outlier fraction and photo-z 
accuracy do not appear to have a redshift dependence 
(for z ~ 0-4). The photo-z quality is comparable to that 
of some recent photo-z's derived for AGNs at brighter 
X-ray fluxes, and is significantly (> 50%) improved 
compared to that of the photo-z's for the X-ray sources 
in the « 1 Ms CDF-S or w 2 Ms CDF-N. 
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TABLE 1 

Summary of the Likelihood-Ratio Matching Parameters and Results 



Catalog 


Del. Thresh 


Depth 


Solid Angle 


% 






R 


C 




Md 


A'NoID 




Atalse 


False % 


Recovery % 


X-O % 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


(18) 


WFIR 


la 


27.3 


1420 


30 345 


0.1" 


0.55 


0.97 


0.73 


462 


344 


118 


2 


32 


9% 


67.5% 


1.0% 


19 


GOODS-S z 


0.6(T 


28.2 


160 


33 955 


0.1" 


0.80 


0.93 


0.82 


311 


259 


52 


8 


23 


9% 


75.9% 


0.7% 


220 


GEMSz 


1.7(7 


27.3 


830 


22 016 


0.1" 


2.35 


0.97 


0.67 


462 


312 


150 


5 


14 


4% 


64.5% 


1.4% 


89 


MUSIC K 


1 CT 


23.8 


140 


13 595 


0.1" 


1.30 


0.93 


0.84 


262 


223 


39 


4 


8 


4% 


82.1% 


1.6% 


14 


MUSYC K 


23.5/arcsec^ 


22.4 


970 


6998 


0.2" 


0.85 


0.99 


0.70 


462 


326 


136 





12 


4% 


68.0% 


4.5% 


9 


SIMPLE 3.6 urn 


2 o- 


23.8 


1640 


22 095 


0.3" 


0.20 


0.99 


0.89 


462 


414 


48 





32 


8% 


82.7% 


1.7% 


12 


VLA 1.4 GHz 


5 a 


19.9 


1170 


338 


0.1" 


2.25 


0.99 


0.20 


462 


94 


368 





0.5 


1% 


20.2% 


27.7% 


83 



Note. — Col. (1): OMR catalog. References.— WFl and GOODS-S: GiavaUsco et al. (2004); GEMS: Caldwell et al. (2008); MUSIC: Grazian et al. (2006); MUSYC: lTk"vlor et alj 
(2009); SIMPLE: Damen et al. (2010); VLA: Miller et al. (2008). Col. (2): The miniinum threshold used for source detection in the catalog. Note that in some cases multiple searches 
have been peii'ormed with different thi'esholds for deblending purposes. Col. (3): Catalog depth in AB magnitude. The AB magnitudes for radio sources were converted from the radio 
flux densities, m(AB) = -2.5 log(/„)— 48.60. Col. (4): Catalog solid-angle coverage in units of arcmin- . The GEMS catalog includes GOODS-S vl data. All the catalogs cover the 
entire Ri 2 Ms CDF-S (ft: 436 arcmin^) except for the GOODS-S and MUSIC surveys. Col. (5): Number of OMR sources in the ft: 2 Ms CDF-S region. Col. (6): 1 u positional error 
of the ONIR sources. Col. (7): Threshold value for the likelihood ratio to discriminate between spurious and real identifications. The thi'eshold value is catalog dependent, and generally 
scales with the typical values of likelihood ratios (see Eq. 1), which usually increase when the catalog depth or positional errors decrease. Col. (8): Sample reliability. See j|2.2l for 
details. Col. (9): Sample completeness. See j|2.2l for details. Col. (10): Total number of X-ray sources that are within the coverage of the ONIR catalog. Col. (11): Number of X-ray 
sources identified with at least one ONIR counteipart in this catalog. Col. (12): Number of X-ray sources not identified, whi ch equals Col. (10) minus Col. (11). Col. (13): Nu mber of 
X-ray sources identified with two ONIR counterparts in this catalog. Col. (14): Expected number of false matches. See j|2.3| for details. Col. (15): False-match probability. See j|2.3| for 
details. Col. (16): Counterpart recovery rate, defined as the expected number of true counterparts (Nm—N^Ax) divided by the number of X-ray sources (Nx)- Col. (17): The fraction of 
ONIR objects that are detected in the X-ray, defined as the expected number of true counterparts (Md "A'palse) divided by the number of ONIR sources in the CDF-S region (A'o). Col. 
(18): Number of primary counterparts selected from this catalog. See i|2.2l for details. 
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TABLE 2 

List of the Primary Counterparts 



ID 




ONIRID 


ONIRRA 


ONIR Dec 


To 


Dis 


Dis/<T 


Rc 


CAT 


MAG 


Multi 


Note 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


1 


0.77 


24484 


52.8924 


-27.8347 


0.1 


0.89 


1.14 


0.88 


GEMS 


24.14 







2 


0.94 


21150 


52.8990 


-27.8597 


0.1 


0.29 


0.31 


0.97 


GEMS 


22.27 







3 


0.65 


24081 


52.9171 


-27.7962 


0.1 


0.44 


0.68 


0.98 


GEMS 


22.26 







4 


0.55 































Edge 



Note. — Table 2 is presented in its entirety in the electronic edition. An abbreviated version of the table 
is shown here for guidance as to its form and content. The full table contains 13 columns of information 
on primary counterparts to the 462 X-ray sources. 



TABLE 3 

List of the Secondary Counterparts 



ID 




ONIR ID 


ONIR RA 


ONIR Dec 


Co 


Dis 


Dis/o- 




CAT 


MAG 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


21 


0.51 


24I5I 


52.9602 


-27.8768 


0.3 


1.38 


2.34 


0.95 


SIMPLE 


21.40 


25 


0.39 


35263 


52.9648 


-27.7648 


0.3 


1.56 


3.16 


0.65 


SIMPLE 


21.99 


31 


0.62 


23319 


52.9752 


-27.8349 


0.1 


0.73 


1.16 


0.63 


GEMS 


21.23 


41 


0.43 


30803 


52.9906 


-27.7022 


0.1 


0.34 


0.78 


0.98 


GEMS 


23.39 



Note. — Table 3 is presented in its entirety in the electronic edition. An abbreviated 
version of the table is shown here for guidance as to its form and content. The full table 
contains 1 1 columns of information on the secondary counterparts of 72 X-ray sources (77 
lines). A high Rc value for the secondary counterpart indicates that the primary and secondary 
counterparts were identified in different catalogs. 



TABLE 4 
Photometric Catalog 



ID 


GALEX 












VIMOS 








NUV 


Err 


Flag 


FUV 


Err 


Flag 


u 


Err 


Flag 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


Median 


23.39 


1.32 


99 


23.87 


1.66 


62 


25.81 


3.79 


440 


1 


-26.89 


99.00 


1 


-27.25 


99.00 


1 


-28.53 


99.00 


1 


2 


23.36 


0.08 


1 


-27.24 


99.00 


1 


99.00 


99.00 





3 


-26.90 


99.00 


1 


-27.23 


99.00 


1 


22.32 


0.01 





4 


99.00 


99.00 





99.00 


99.00 





99.00 


99.00 






Note. — Table 4 is presented in its entirety in the electronic edition. An 
abbreviated version of the table is shown here for guidance as to its form and 
content. The full table contains 181 columns of photometric data for the 462 
X-ray sources (Rows 2-A63, including 16 sources with no counterparts). The 
first row includes the median magnitude information for each band. All the pho- 
tometric data are given in units of AB magnitudes, and are set to "99" if not 
available; upper limits are indicated by the "-" sign. 



TABLE 5 
Photometric Redshift Catalog 



ID 


Spec-z 


Photo-z 


zLower 


zUpper 


#Detect 


#Filter 




PhOtO-Z;||| 




Op. Det. 


zRef 


(!) 


(2) 




(4) 




(6) 


(7) 


(S) 




(10) 


(11) 


(12) 


1 


-1.000 


1.102 


1.045 


1.114 


13 


16 


1.8 


-1.0 


-1.0 


1 


-1 


2 


-1.000 


1.626 


1.620 


1.633 


26 


27 


23.6 


-1.0 


-1.0 


1 


-1 


3 


2.709 


2.714 


2.712 


2.715 


26 


28 


8.3 


-1.0 


-1.0 


1 


8 


4 


-1.000 


-1.000 


-1.000 


-1.000 








-1.0 


-1.0 


-1.0 





-1 



Note. — Table 5 is presented in its entirety in the electronic edition. An abbreviated version of the table 
is shown here for guidance as to its form and content. The full table contains 12 columns of information 
for the 462 X-ray sources (including 16 sources with no counterparts and six stars). 



